VIBRATIONAL PROPERTIES OF DISORDERED 

BINARY ALLOYS 


By 

RAVI PRATAP SINGH 




' DEPARTMENT OF PHYSICS 

INDIAN INSTITUTE OF TECHNOLOGY KANPOK 

NOVEMBER, 1982 





VIBRATIONAL PROPERTIES OF DISORDERED 

BINARY ALLOYS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 



By 

RAVI PRATAP SINGH 


CO 

department of physics 

INDIAN INSTITUTE OF TECHNOLOGY KANPUB 

NOVEMBER, 1982 



1 - lUt' 1984 


fISli 


PHV- /as? ~D- Siw- 



iSrSTTT 


rT" 





11 . 


CERTIFICATE 

This is to certity that the work in this thesis 
entitled: "VIBRATIONAL properties of disordered binary ALLOYS 
has been carried out by RAVI PRATAP SINGH under my supervision 
NO part of this work has been submitted elsewhere for a degree 


November/ 1982 


Abhi^t Mookg^ee 
Assistant Proressor 
Departnent o£ Physics 

Indian Institute of Technology 




POSr GRADUA TE OFFICE 
This thi*''vls hrns bc.tm apptovcti 
for thr .-wanl of the Degw of 
Doctor of Philosophy (Ph.D.) 
i in actofu.tncc with the 
rcgulatioiis ol Uic iiuiian 
Jnsiitutc of Technology ■ anpuf 

1 l l /s/Jy ^ 



iii 


ACKNOWLEDGEMENTS 

The able guidance, constant help and wise counsel of 
Dr. Ahhijit Mookerjee is gratefully acknowledged. His patient 
care made it possible for me to carry out and complete this 
investigation. At times I could afford to falter for I knew that 
he was around to help me and he did it. It is embarrasinyly 
impossible to list out what all I got from him. Working with him 
has been an unforgettable experience of life. 

My sincere thanks are due to Dr. M. Yussouff for many 
helpful discussions at various stages of this work. His constant ‘ 

I 

encouragement and appreciation about the progress of my work | 

helped me in achieving the goal without hurdles.. I wish to thank j 

[ 

Dr. M.J. Kelly of University of Cambridge for helpful discussions j 

1 

and Dr. R. Haydock of University of Cambridge for providing the | 
xaseful programmes. I am thankful to Professor D.C. Khan for 
talcing interest in the progress of my research work. 

I thank my friend Mr. Pankaj Joshi for his critical 
cojiroents from time to time. 

f 

f 

I must thank my friends Dr. A.K. Nigam,- J.K. Sharraa, 

V.K, Srivastava, Vasundhra Choudhxry and Messrs A.K. Upadhyaya, 

Arun Agrawal, B.P. Singh, D.S. Misra, V.D. Chafekar and many 
others who helped me in more than one ways during the course of 
* this work.. 

It is a pleasure to thank Dr. M. Misra for his timely ; 

help in many ways at the time of completion of this thesis. ; 



My thanks are also due to the Director (Personnel)/ Oil 
and Natural Gas Commission/ Mr. M.M. Goswami (Supdg. Geophysi- 
cist (s)) and Mr. M.B. Singh/ Sr. Geoph;ysiGist (s ) / ONGC/ Nazira 
for providing me leave for completing this work. 

I am grateful to Council of Scientific and Industrial 
Research/ New Delhi for financial support in the form of Junior 
Research Fellowship for 2 years during this research work. 

I thank Mr. R.N. Srivastava for neat typing, Mr. J.S. 
Sharma for tracing the figures and Mr. H.K. Panda and Mr. L.S. 
Rathore for cyclostyling this thesis . 

I take this opportxinity to pay homage to my dear friend 
late C.P.K. Reddy whose untimely death has created a void in my 
life . 

My warmest thanks go to my dear friend Gyanesh and to my 
younger brother Rajive whose company made a homely atmosphere for 
me during my stay at I IT Kanpur, ‘ 

In the end/ I would like to appreciate the importance of j 
computer DEC-10 of IIT Kanp\ar without which it was q-uite difficul-^ 
to achieve my goal, I 


(Ravi Pratap Singhj 



CONTENTS 


Chcfpter 

LIST OF TABLES 
LIST OF FIGURES 
SYNOPSIS 

I . Introduction 

1^1 Lattice Dynamics of Disordered System 

1«2 The Harmonic Model 

1-3 Lattice Green ' s Functions 

1.3.1 Equation of Motion of the Green's 
Function 

1.4 Lattice Symmetries 

1.4.1 Reduction of the Force Constant 
Matrix for Cubic Symmetry 

1.5 Green's Function and Dispersion Curves 

1.5.1 Born -von Karman Models in Cubic 
Lattices 

1.6 Experimental Quantities 

1.6.1 Density of States 

1.6.2 Inelastic Neutron Scattering 

II, The Coherent Potential Approximations 

2.1 Introduction 

2.2 Virtual Crystal Approximation (VCA) 

2.3 The Average T-matrix Approximation (ATA) 

2.4 Single Site Mass Coherent Potential 
Approximation (ICPA or MCPA) 

2.4,1 Failures of MCPA 

2.5 Generalisation of CPA 

2.6 Cluster Embedding 

2.7 'Augmented Space* Formalism 

2,8. Lattice Green's Function in 'Augmented Space' 

2.9 The Graphical Method 

2.10 CCPA Self-consistent Medium 

2.11 Spectral Functions and Density of States 

2.12 Results on Linear Chains 


Page No 

vii 

viii 

xii 

1 

1 

3 

5 

7 

9 

13 

20 

24 

25 

25 

26 

28 

28 

28 

29 

31 

34 

35 
37 
39 
45 
48 
50b 
59 
62 



III. Application of CCPA to FCC Structured Transition 


Metal 

Alloys 


3.1 

Introduction 

69 


3-2.1 Force Constant Matrices 

71 


3.2.2 Lattice Green ' s Functions 

72 

3.3 

Ni-Pd Alloy 

79 

3.4 

Ni-Cr Alloy 

87 

3.5 

Ni~Pt Alloys 

92 


3.5.1 Nig^Pt^^ 

92 


3.5.2 Ni^^gPt^Q 

98 


3.5.3 

104 

3.6 

Conclusion 

106 

IV- Concluding Remarks 

107 

References 

111 

Appendix A 

Recursion Method 

A-1 

Appendix B 

Dielectric Function, Static Screening, and 
Non~Linear Response 

B-1 



LIST OF TABLES 


Table 

3.1 

3.2 


Page No. 


General Properties of FCC Structured Metals 

(Ni, Pt, Pd and Cr) 70 

Recursion-Coefficients of FCC (Ni) Lattice 70 



Figure 

1 .1 (a-ti ) 

(a) 

(b) 

1 . 1 ( c-e ) 

(c) 

(d) 

(e) 

1 . 2 ( a-c ) 
l,3(a--d) 

(a ) 

(b) 

(c) 

(d) 


1.3(e) 

1,4 (a) 
(b-G ) 
(d-f) 


LIST OP FIGURES 


Page No 

Point symmetries of Bravais lattices 12 

One -dimensional lattice: unity and inversion, 

X = + X; inversion indicated by arrows 

Two-dimensional square lattice : all 8 opera- 
tions which leave the square invariant 
(X, Z) = (X^, X^) = (+ X^, ± ^3^' 

Cubic symrnetry operations ' 12 

Fould fold axis j| cubic edge, '[^lOo] 

Three fold axis | ] body diagonal Llllll 
Two fold axis ] | face diagonal (^lioj 

Coupling in j^OOj direction 14 

Coupling matrices for cubic symmetry 16 

R“ = Of self-coupling 0^°^ is isotropic (scalar) 

^ I![l00], (1.33), contains 2 para- 

meters, one longitudinal spring f^^ and one 
isotropic transversal spring f^. 

R- II [ill], (1.34), as in (b), 

f^=a+23, ft=®~P 

R^ II [ 1103 , (1.35), contains three 

parameters , one longitudinal ' spring ' t and 
two different transversal springs f ^ 
fl = c(+6,ft, = a-B,f^ = Y 

Nearest neighbour model for a FCC lattice with 16 

three parameters a, 6 / y or f^^, f^,/ 

A sxibstitutional defect representation 19 

Transversal springs coupling representation 

Violation of rotational invariance for a 
defect model 





1.5 


2.1 


2.2a 


2.2b 

2.2c 


2.3 


2.4 


2.5 


2.6 


3.1 


3.2 


3.3 


Dispersion curves of FCC Ni metal from a 
Born->von Karman model fit with nearest neighbour 
springs in the main symmetry directions and the 
exx^erimental points from Birgeneau et al. (1964) 

The augmented space octagon decorating a bond 

Renormalization of the bonds by augmented 
s pace octagons 

Renorn'alization of the ocuagon by the medium 

The augmented space 24-gon for a 3CPA calcu™ 

1 at ion 

Linear chain phonon density of states for small 
mass and force constant disorders 

Comparison between the embedding in a ICPA 
iTiedium self-consistent Tsukada cluster medium 
for the linear chain 

50-50 alloy with a light component with 
enhanced force constants. Comparison with the 
exact results from computer calculations 

A dilute/ light impurity with enhAnced force 
constants showing an impurity band with consi- 
derable structure 

j Dispersion curve's for FCC Pt metal from a 
Bom— von Karman model fit with nearest neighbour 
springs in the main symmetry directions and the 
experimental points from Dutton et al. (1972) 

Dispersion curves for FCC Pd metal from a 
Born-von Karman model fit with nearest neighboxir 
springs in the main symmetry dixectiona and the 
experimental points from Muller and Brockhouse 
(1972) 

Dispersion curves for hypothetical FCC Cr 
metal from BCC elastic constants. Here maximum 


23 


51 

53 

53 

53 

63 

64 


66 


67 


73 


74 


75 



3.4 


frequencies in all symmetry directions are 
matched with BCC values (Feldmann, 1970) 


77 


3.5. 


3.6a 


3.6b 


3.7 


3 . S 


3.9 


Frequency distribution for pure Ni metal from 
the recursion method. The frequencies on the 
abscissas are in THz and the units on the 
ordinates are chosen so that the distribution 
functions ore normalized to unity 

Frequency distribution for Ni metal from 
oxpirimontal data (K-space. method) . Hero also 
the ordinate is chosen so that the distribution 
function is normalized to unity 


Spectral functions for Ni^^Pd^^ in 2CPA in high 

symmetry directions. The frequencies on the 

abscissas are in THz end the units on the 

ordinates are chosen so that the functions are 
/ 

normalized at the maximum peak value 


Spectral functions in 2CPA for 55^^45 
direction. The frequencies on the 
abscissas arc in THz and the units on the 
ordinates aro chosen so that the functions are 
normalized at the maximum peak value 


The width functions for ^.i^^Pd^^ in MCPA and 
2CPA. The experimental points are taten from 
Kamitakahara .and Brockhouse (1974) 


Dispersion cxrrves for Ni^^Pd^^ in MCPA and 
2GPA in high symmetry directions. The full 
curves are for MCPA; broken curves are for 
2CPA 

Frequency distributicai for Ni^j^Pd^^ in MCPA 
and 2CPA. The ordinates are chosen so that 
the distribution functions are normalized to 
unity. The full curve is for MCPA; broken 
curve for 2CPA 


78 


80 


81 


83 


85 


86 



3.10 


3.11 


3.12 


3.13 


3.14 


3.15 


3,16 


Some spectral functions for 2 

MCPA and 2CFA. The ordinates are normalized 
at the maximum peak value. The full ciirves 
correspona to MCPA; broken curves are for 2CPA 


Dispersion curves for Wi-^Cr, ^ HCPa and 
2CPA in high symrrietry directions. The full 
carves are for MCPA; broken curves are for 
2CPA 


88 


89 


Frequency distribution for NiygCr ^2 liOPA 
and 2CPA. The obdinates are chosen so that 
the distribution functions are normalized to 
unity. The full curve is for MCF/\; broken 
ciirves are for 2CPA 


91 


Dispersion curves for ^^-^ 95^^05 MCPA and 
2 CPA in [ 00 direction. The full curves are 


for MCPA and broken curves arc for 2 CPA. 
circled positions are enlarged 


The 


Spectral functions (lino shapes) for Ni^jgPtQ^ 
in 00 C direction for MCPA and 2CPA near 
resonance mode. The ordinates are normalized 
at the maximum peak value 


Disjxrsion curves for Mig^PtQ^ in MCPA and 
2 CPA in [zcoj and directions. The full 

curves arc for MCPA; broken curves are for 
2CPA 


Frequency distribution for Ni^gPtQg in MCPA 
2CPii. The full exorvo is for MCP/i; broken 
curve for 2Cr-A- The ordinates ore chosen so 
that the distribution functions are normalized 
to unity. In the inset are given the experi- 
mental density of states by Tsunoda et al, with 
the partial density of states (broken ciorve) 
of Ni in the Nig^PtQ^ calculated with the MCPvk 


93 


94 


95 


97 



3.17 


3.18 


3.19 


3.20 


xii 


Dispersion curves for Ni^gPt^Q in MCPiv and 99 

2CP7i. along the- high symmetry directions. The 
full curves are for broken curves are 

for 2C17i 

Frequency distribution for Ni^Qpt^Q in MCFii. 101 

and 2Ci:7i. The ordinates are chosen so that the 
functions ore normalized to unity. The full 
curve is for MCPA and broken curve for 2CPA. 

In the inset are given the expcsrimental density 
of states by Tsunoda et al, with the partial ! 
density of states (broken curve) of Ni in the 
Ni- 7 QPt 3 Q calculated with the MCPA ' 

Dispersion curves tor hi^Qit^Q in hCPiv ‘ana ‘ 103 

2CFA in high symmetry directions. The full 
curves are for MCIA; broken curves are for 

2cr;i 

Frequency distribution for Ni^Qpt^Q in MCP^* 105 

and 2CI.V. The ordinates are chosen so that 
the functions are normalized to ■unity. The 
full exorve is for MCIFx and broken curve for 
2CI7V. In the inset -'‘xe given the experimental 
density of states by Tsunoda et al. with the 
partial density of states (broken cuirve) of 
Ni in the ^iso^^SO calculated with the MCPA. 



xill 


SYNOPSIS 


When atoms of different masses (site energies) and diff- 
erence interatomic force constants (overlap-integrals) alloy 
substitutionally in a crystal^ there are/ in general three types 
of changes in the system — Hamiltonian diagonal 0 . . , off- 
diagonal (where j is neighbouring site of 1) and onvir- 

onmental Unlike electronic systems / such types of disorder 

J-Ki 

are coupled in a vibrational system by the general translational 
and/or rotational symmetry of the system (Chapter I ) ^ rendering 
its study extremely cumbersome . 

The first successful theoretical study of substitutionally 
disordered vibrational system was due to Taylor (1967) known as J 
mass-defect coherent potential approximation (MCPA or ICPA). The 
MCPA is a self-consistent single-site approximaticxi for alloys ; 
with diagonal disorder. The need to go beyond MCPA has been felt | 
for a long time. Several suggestions have been put forward about | 
the ways of generalizing the MCPA and ways of incorporating off- 
diagonal disorder, environmental disorder and short range-order 
(SRO) (Briefly reviewed in Chapter II of the thesis). CalcxxLa- | 
tions based on the^e suggestions suffer from two main drawbacks. 
Firstly, for any exactly solvable realistic irKxiel, in which all 
types of disorder are coupled, the corresponding algebraic equa- 
tions become tediously unwieldy, particularly on a three-dimen- 

I 

sional lattice. Secondly, any tractable or further simplified 
approximation always leads to the unphysical averaged Green *s 
function in the strong scattering regime. The essential Herglotz 
property of the Gref's fxonction is violated. 
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The "augmented space" of MooJserjee (1973) preserves the 
correct rierglotz x^roperty of Green's function. This forroalism is 
discussed in detail in Chapter II. The present work which is 
based on " augmented space" offers an unambiguous tractable / self- 
consistent clust-^r gon-..-ralieation (CCPA) of MCPA, which retains 
the Herglotz prop-.5rty at all frequencies for all ranges of dis- 
order. The general formalism of the CCPA is also discussed in 
Chapter II . A graphical me-choo coupled with the recursion tech- 
nique of Haydock (1972) enables us to work on a realistic three- 
dimensional lattice , 

The most important contribution of this work is the . 

generation of a self-consistent self -energy- This produces the | 
vibrational properties (viz. dispersion curves, density of states, 

I 

f 

neutron scattering cross-section, etc. ) more accurately than that j 
produced by MCPA- The results of self-consistent cluster calcu- 
lations for the density of states for linear chain are reported 
in Chapter II. The band cages are much better reproduced, as well j 
as the struct-ure near the bond edges, especially in those cases 
where there is an impurity bond with consider cible stiructure. 

In Chapter ill our cluster - CPA formalism is applied for 
the calculation of vibrational properties of random transition 
metal alloys c.g. Ni^_^Pt^, Ni^__^Pd^, Ni^__^Cr^ (0 < c < 1) where 
effects of off-diagonal disorder and clusters are expected to be 
important. These alloys form a continuous series of PCC solid 
solutions and thdir experimental resxxlts are available to compare 
with the theoretical ones. Bosi et al. (1980; 1978) have applied 
^!CPA to these alloys. They have treated the diagonal disorder 
in MCPA but off-diagonal disorder in virtual crystal approximation 
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(VGA). This is not a consistent procediore. The fit with experi- 
ment is not satisfactory and clustering has been forwarded as 
being responsible. We have taken off-diagonal disorder into 
account and the cluster coherent potential medium is made self- 
consistent. Our results ar^; closer to experimental results. The 
trend of changes due to clusters is in the right direction. 

Last but not the least is the work related with the 
screening density around an impurity of unit charge given in the 
Appendix (R.P. Singh and M. Youssoxrff, Phys. Stat. Sol, (b) 101 ^ 
K93 (1980)). 

In the concluding Chapter IV, the possible extension and 
future usefulness of the cluster CPA have been discussed. 


‘ii 



Chapter I 


Introduction 

1.1 Lattice Dynamics of Disordered Systems ; 

The study of matter with ri^dom composition nas for many 
years received much less attention from physicists than better 
definable crystalline materials. Nevertheless, some important 
iresults on average permittivity of a heterogeneous medium were 
obtained very early in certain disordered systems (Rayleigh, 1892). 
Apparently the first work devoted to a calculation of the frequ- 
ency spectrum of a disordered three-dimensional crystal was that 
of Lifshitz and Stejaiova (1956), These authors obtained the 
leading terms in the expansion of the spactrim in powers of the 
concentrdtioh of the minority species and determined the' shift of 
the maximum frequency of the unperturbed crystal due to the 

introduction of a finite concentration of impurities. Their I 

i 

results were not applied to the nxamerical calculation of the I 

spectrum of any crystal model, • j 

Aii essentially analogous method for expanding ■the spec- 
trum in ptowers of the impurity concentration, based bn the use 
of Green *s fxinction was proposed for the linear chains by Davies 
and Danger (1963), This method was subsequently generalized to 
the three-dimensional disordered crystals and the calculabions 
•were made seif -consistent by Taylor (1967) ensuring that his 
approximations to the frequency spectrum should be independent of 
which type of atari is regarded as host and which is regarded as 
the impurity. One of the consequences of ensuring ■that this 
svmmetrv is present in the theory is that the impurity bands which 



are predicted do not overlap the maximum frequency of the monoa- 
tomic light -atom lattice. Within this formalism - widely known 
as Mass-defect Coherent Potential Approximation (MCPA) - however/ 
the randomness in the force constants could not be taken into 
account. Moreover being a one-site approximation/ effects of 
random clustering also could not be trck-n into account. 

Despite these approximations/ the MCPA theory has baen 
useful in explaining the experimental results of the alloys with 
small mass-aisoraer such as Nb-Mo, Cu-Zn/ Cu-Ni (Grionew-ald/ 1980)v 
In case of large mass -disordered alloys such os Pd-Pt/ I'lb-Ta, 

Ge-Si etc./ the largo resonances and localized modes observed 
experimentally were found in general to be in fair agreement with 
MCPA provided that the disorder was essentially in the masses 
only and not in the bonding (Bosi et al, 1978). The agreement 
was \insatisf actory for the alloys such as Cr-W/ Cu-Au, Ni-Pd/ , 

Ni— pt. The expcriiuants have shown both resonant and localized 
modes. For resonant modes it was found that the frequency shifts 
were branch dependent for Cr-W and Cu-AU/ contrary to the expec- ■ 
tat ion of MCPA theory (Cunningham et al/ 1970; Svens son and 
Kamltakahara/ 1971). It was foxind that the zero frequency shift 
is at higher energy, than the resonance in Cr-W and Cu— Au alloys. 
For Ni-Pd tho phonon width dependence on frequency is different 
for various branones and in disagreement with MCPA calculations 
(Kamltakahara and Brockhouse/ 1974), In the case of Ni-Pt, the 
energy splitting of the phonon dispersion curves at the resor- 
nance mode frequency has a branch dependence in contrast to the 
MCPA prediction. Moreover, the experimental phonon density of 
states deviates from the MCPA calciilations for the samj^ies with 



concentration C 0,30 cind the difference is most conspicuous 
for the maximum jehonon frequency which decreases faster with incre- 
asing Pt concentration than the theory predicts (Tsionoda et al, 
1979). 

Such deviations urge the need of incorporating other 
disorders viz. off-diagonal and environmental, which are coupled 
in the case of phonons, by the genc.ral translcitional and rota- 
tional symmetry of the system. In the following work we shall 
concentrate on substitutionally disordered binary alloys with 
off-diagonal disorder. In view of this, we have first reviewed 
the requirements of detailed study of disordered alloys in this 
chapter. In Chapter II, we review the earlier methods (i.e. 

MCPA, etc. ) with their drawbacks and discuss possible lines of 
generalisation with the self consistent cluster coherent potential 
approximation (CCPA) based on "augmented space" of Mookerjee 
(1973). This formalism (CCPA) has been applied to transition 
metal alloys e.g. Nij^__^Pt^, Ni^__^Cr^ and Ni^__^Pd^ which fall in 
the category of substitutionally disordered alloys in which the 
effects of off-diagonal disorder <:ind clustering are expected to 
be important (Chapter III ) . In the concluding Chapter IV, we 
have discussed the main achievements of CCPA formalism in the 
stxady of vibrational properties cind its various possible appii- 
"cations to similar problems. 

1.2 Iha Harmonic Model 

In the harmonic approximation, the Hamiltonian for the | 

j 

vibrational motion of on ordered crystal is - j 

P^(l) 

^ + I . d/l' ^ (l.l) 

la 2M° 1 a “ “ f 

^ I Cli f 



with tho force constant matrix ^ being translation ally 
invariant. Here Ua(l) and P„(l) the Cartesian coordinates of 

the displacement and momentum operators, respectively, of an atom 
in the unit cell at R(l), coordinate, where Mq is the mass of 
that atom. We t-^ke the same form for the alloy i.e. 

A. r ^ 

^ ^ airTT) ^ u^^(i).j2>„c, . (i,i'),u^, (1* ) (1.2) 

1C£ Ot 1 Cl 

1' O' 

Note that the mass now depones upon 1 and that 0 is not 
translationally invariant. The Heisenberg equations of motion 
for these operators are 

ih'—xiad) = t u^(l,t), h 3 = M^^(l) p^(l,t) (1.3a) 

ih 1^ P^d) = [;p^(l,t), hJ = - Zi2fca,dd‘)u„,d',t) (1.3b) 

So that 

M^jd) u (t) = ~ J (1,1-) u (l',t) (1.4) 

?t^ 1' a' ““ 

When Potirier transformed to frequency w, this expression 

beccxmis 

Z Cm„( 1) S « (1,1') - J3 u„,Cl') = 0 (l.S) 

Ql 2 t 

In the perfect crystal, the characteristic frequencies 
and normal modes may be recidily obtained because of the transla- 
tional symmetry (Maradudin et al, 1971). But in the disordered 
crystal the translational symmetry is lost. Nevertheless, the 
equation of motion (l.S) in the disordered crystal can conven- 
iently be expressed in terms of that for the perfect crystal where 
it becomes 



0 


I 


aa 


11 


^aa ’ 


(1.1‘) - V 


aa 


( 1 , 1 -)] 


u 


( 1 ) 


( 1 . 0 ) 


where 




*^0 ■'«! 


■ '^11> 




with 


Mad) = M a - 



(1 -£,^ (D) 


(1.8a) 


and 


Paa.(ld') 



(1,1' ) + A ^aa' ^ 


(1.8b) 


This problem Is formerly identical to the problem of 
non-interacting electrons, in disordered crystals (Elliott et al, 
1974). Within the harmonic approximation, the mow ledge of fr-g- 
uency and spatial extent of each normal mode is necessary, which 
could bo obtained for the translationally invariant crystal by 
the use of symnetry transformation. In on Imjerfeot crystal 
this procoduro is no longer of use and other means of attach on 
the equation of motion (1.5) are required. It is a fact that 
the experiments on disordered crystals (viz. neutron scattering, 
optical absorption, ate.) simply measure a correlation or response 
function rather than the properties of particular modes (Van 
1954). A well defined formaUsm which leads Immediately to these 
quantities without the necessity to solve for the normal modes 
is that of thermodynamic Green's functions which we use and 

discuss below. 


1.3 _^ttice Green's Functions 

The required function of the lattice dynamics is almost 

always the displacement-displacement correlation function 



C30 


iwt 


S,,a.(l/l';w) = / <u^(l,t) . u^,(l'.0)>^ e— caw 

- 00 


(1.9) 


where <... cienotes a thermal average. 

Such correlation functions are readily calculated via Green's 
functions (Fetter and Walscka# 1971). We choose to use Green's 
functions given by 

Gctai(lyl';w) = (^) j u^,(l',0)» e’^'^^dt (1.10a) 

- 00 


where 

«A(t); B » = tie (+t) <[A(t), Dj>^ 


(1.10b) 


is a double time Green's function (Zubarev, 1960). The upper 
and lower signs give respectively, the retarded and advanced 
Green's fxjinctions with 0(t) being the usual Heavisic3e unit step 
function. 


Thfjsc Green's functions are related to the corirelation 
funcrtion (1.9) via the "spectral representation" 


G 


aai 


(l,l';Z) = (:^) / 


CD 

- CD 


(1 - Saa. (l.l';w) 

Z -- W + i*e 


dw 


( 1 . 11 ) 


where 3 a -h/iCgT, A more explicit relation between these fxoncticjns 
is obtained" from the discontinuity of the Green's fxonction across 
the real axis i.e. 


G(w + io^) - G(w - io"*") 


2i Im G(w) 


(^) (1 - S(w) 


( 1 . 12)1 


Prom the analytic form (1.11) it is clear that G(w) is 
analytic in the upper (lower) half complex 2 plane and thus the 
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real and imriginary parts of G(w) are related by Kramers -Kronig 
relations 


Re G (w ) 



Im G (w » ) 
w v; * 


dw ‘ 


(1.13a) 


and 


-j CD Re G ( w * ) 

Bn G(w) = (- -i) p / -- — dw' 

TT w - w 


(1.13b) 


where P / is the principal part of the integral. Since displa- 
cement operator UQ(l;t) is oven under time reversal, we find 

G(-w+i<5)=:G(w+i<3) (1.14) 


So that the jriilbert transform (1.12) may be folded into integrals 
over only the positive frequencies 


Re G (w ) 


o 


00 


v; ' Jm G (.w ' ) 


dw ' 


(1.15a) 


w 


w ' 


and 

Im G(w) = (- ^) 

The major experimental quantities (viz. phonon disper- 
sion curves, density of states, etc.) can be determined if the 
Green's functions can be found from the equation of motion whidi 
we discuss next. 

1,3,1 Equation of motion of the Green's fionction 

Differentiating the formal Green's function (1.10) with 
respect to time gives the equation of motion 


00 Re G(w' ) 

/ o dw' 


w 


w 


.2 


(1.15b)! 
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±R (■—) « U,. (l';0)» 


= 2fr>^ (t) < [u^(l), U^. (1* )J > 

+ «[u^(l;t)^ Hj; u^, (l';0)» 

= 0 + ±K « p^(l,t)/My(l); u^,,(l',0)» (1.16) 

This contains a different Green's function/ so proceeding 
to the next equation 

ih (~) «p^(l/t)/M^(l); u^,(l';0)>> 


2 ^^ (t) <[p (1)/M (1); u 

a ct cx ^ 

+ «[p^Cl;t)/M^(l)/ hJ; u^.(1‘;0)» 


2 IT ifi X fi. \ A 6 


m„(iT 


^an ' 11' 


c-y - V T' ^ jZ „(1/1") . «ua„(l"/t); u„, (l';0)» 

(1.17) 


Defining 




« u^(l;t)/ u^, (l‘;0) » 


(1.18) 


We find 


2 ^rta' ^ 

gt 


(e hi. h. - 4„ 


(1.19) 


The Fourier transform to frequency gives 



(1,1') + Z ,, 

a'>i" 


( 1 . 20 ) 




M^(l) w (l,l';w) = 


G „ , (1" /!' 7W) 
a” a' ' ' 


Written in the full matrix notation 


M G(w) = r + ]^ G{w) 
or 

G(w) = [Mw^ - (1.21) 

The Green's function for a perfect crystal described by 
the Hamiltonian (1.1) satisfies a similar equation 

G°(w) = (1.22) 

Hence on introducing the perturbation matrix 
y(w) = (M° - M) +[^ ~ (1.23) I 

We can obtain the basic Dyson equation relating G(w) 
and G° (w ) , 


G(w) = G°(w) -t- G°(w) V(w) G(w) (1.24) ; 

mmm mm* mmm mm, ^ f 

f 

f 

I 

Solution to this equation will be described in Chapter II. * 
However it is clear that a detailed knowledge of G°(w) is required ; 

and we, hence/ direct our attention to this problem. We first ; 

1 

. O , V ! 

discuss lattice symmetries and their utility in evaluating G (w) ! 

i 

and G(w) as well. i 

[ 

1.4 Lattice Symn^tries : 

Synmetry operaticxiS/ S/ consist of a special translation | 

I 

and/or rotation s(T/r)/ which leave not only the potential but ! 
also the lattice itself .invariant: [ 



(1.26) 


m m m * 

SR =T + rR = R** 


where rv^ m are lattice site indices, 

o 

m 

i.e. with R / R ^ also covers all lattice sites. These symraetries 
are very important; they allow one to reaucs the number of indep- 
endent force constants and to obtain the eigenvectors of the 
coupling matrix. Therefore we discuss them in detail for three- 
dimensional (3-D) Bravais LEittices. Hero the .translation-rota- 
tional invariance 


Si Si 

j2(( i + ^ + 


m m 

j 2 !( I + rX 4-rS) 


m. 


m 


= j2f( X + r S ) 


(1.26) 


m 


means that the expansion of internal potential U about X with 

S -s 

displacement S is identical to the expansion about X with 

m 

displacements r S , 


m mn n 


n 


2U 


_ _ m Sio^c; _ 

(s . S ) = (r s"", P ^ r S ) 


m -4 ^ 

(s", r” 0 r s'~' 


(1.27) 


which leads to 


mn 

0"" 

mn 

.mm mm 


1 JSs^s 


r 0 


r ; 


r (T r 
i'i ^i'lc’ k'k 


Translational symmetry, m=m+h, r=l means 


mn 


m + h, n 4- h 


(m - n) 


= jZf = iZf 

Inversion about a lattice point, r = -1, m^ 


(1.28) 


(1.29a)| 


- m, xs 


also a symmetry operation for eveiry Bravais lattice. 


r 



(1.29b) 


mn 

iZf— 


+ (m ~ n) 
jZf ~ ” 


mn ran 

^ik ~ ■^ki 


nm 

r‘ ; 


Besides inversion the point symmetries contain rotations 


about lattice points which depend on the type, of lattice. Let 
us discuss (1.28) 





h ho 

r R = r R ^ 


(1.30) 


•otation ^ , 

-.1 h 

r ) r (1.30a) 

h 

consider rotations only which leave R invariant 
h 

of R , then 

-1 

r 0(R ) r 

(1.30b) 

-.1 

r 0 r 

(h) 

form of 0 to make it compatible with the parti- 

Ch) ! 

cxilar point symmetry. After having fotind a compatible 0 , one j 

can find the coupling matrix for equivalent h' from (1.30), I 

I 

Here we discuss only structures of highest syiranetry. I 

Figure (1.1a) shows the linear lattice, with the obvious point | 

symmetries imity and inversion, corresponding tox=(rx)=+X. | 

1 

The symmetries of the plane square, lattice. Figure 1.1b, consist 
of all rotations which leave the elementary square (the Wigner- 

I 

Seitz-Cell) invariant. These are 8 operations (+ X, + Z), 

21, + X). The highest symmetry in three-dirrensional lattices 
is cubic. It contains all 48 rotations (Figures l,lc-e) about 


or in another 
h 

j2((R“") = 

If WG 
for any value 
h 

i^(R“ ) = 

or 

(h) 

restricts the 








the centet of a cube which transf ormjkthat cube into itself : 
(+ X, + Y, + Z) -and pennutations of Z. 


1.4,1 Reduction of the force constant Matrix for Cubic Symmetry 
Let us first discuss = p(o) whichy obviously, is 

invariant against all 48 operations: tne off-diagonal elements 
must vanish, the diagonal elements arc ^qual and is scalar,^ 




■f . 

^ 'ik ' 


1 independent parameter f 


h 


If R ! ! ^lOOj thu rotations (X, + Y, + Z) and (X, 
h 

not change R . The coupling must have the form 


(1.31) 
Z, 4- Y) do 


jZf(R 


h 


[100| ) 




100 


f 0 


0 

0 


0 

0 


"•t 

0 f. 


(1.32) 


2 independent parametersL 

It can be represented by a longitudinal (spiral ) spring 
an isotropic transverse (leaf) spring f^. Figures (1 . 2a-c) , 


^3 


, 100 


,iioo1 


+ f^(l 


% 

\ 

't 




[106] 


« f. ^ P 


1100 ] 


h 


■1 U t ‘t 

with the longitudinal projector in direction R , 

10 0 


(1.33a 


100 


t R > 


0 0 0 
0 0 0 


(1.33]3 


By 


the rotatic^ (- Y, X, Z), where [lOO] [^lOj , one 


obtains from ( 1 . 30a ) 


[ 
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mi 
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0 



a(W) 


Fore# £ oo atoii 0 du« to aun, Arbitrary di* placement « o£ an 
•tow m% Jg» (/ rilO], hmtm ^ » iX, 0, 0). ~ 


/ 

* ✓ 

/ 

pr: ’ r:vj : : 0<t>AAA^WWAWS •» 

b)l Th» iongit'Jdliial cowpoiMwit of loacia to a force fi«i 

in [ lOOJ-dirwction* vfiich emm be represented by e longl* 
tudlnel (spiral} spring fj^* 


* 



5 | 

„.A. 


d) TTe transversal cosipooent of f, •►, leads to a force ft»t 
perpendlculsr to the LlOtg -direction « which can be repre- 
‘'sefited by en Isotropic transversal (leaf) spring f^. 


FIG.H Coupling in[l00] direction 
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0 

fx 


=. f p toioj + f p . 

*1^1 * 


3 1010] 
1 


I 0 0 0 

I 0 1 0 

0 0 0 


this shows how to extend (1.33a) to other <100> couplings, 

h h 


If R 


(h) 


(^lllj all permutations of (X, Y, Z) leave R 


unchanged and 0 has a representation equivalent to that of 
(1.30a) 


- 12 


[m; ^ f £ M _ liiij, [ill] „ 1 


111 
111 
111 
(1.34) 

There are then 2 independent parameters, although the form of 
■(h) 

0 is changed (Figures 1.3a-d). 
h , , _ 

If R Il[ll0l all operations (X, Y, +Z), (Y, X, tZ) are 
admitted. The coupling is represented by 


0 


[IIOJ ^ £^P^[1103 ^ £^^pjll0j ^ ^^pjlio] 


(1.35) i 


t' t' t t 

3 independent parameters 

i.e. a longitiadinal spring and two different transversal springs; 
projects onto a cubic edge 


looil 



" 1 

1 

0 


' 0 

0 

o"" 

„’Liial 1 

^1 ■”2 

1 

1 

0 

! 

piliol = 

' ^t 

0 

0 
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^ 0 
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0- 
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J- # 


,Ciio] 

f 



1 0 
1 0 
0 0 



(1.35a) 


The coupling is still longitudinal and transversal in character, 
but the transversal behaviour is no longer isotropic (Pigxires 
1 . 3a— d ) . 


Now we are in position to coiistruct a force constant model 
for cubic lattices. Because of translational invariance 


I 

h 


(h) 


0 


0 


(o) 


- 1 0 


(h) 


(1.36) 


In cubic crystals, where 0 
Co) „ « tr ^ 


(o) 


is scalar. 


0 


ik 


ik 


I tr ( 0 } 

^ h?^o 


(o) 


ik 


(1.36a> 


r(o) 


A very crude model, where only 0 enters, is the 
Einstein model in which each atari moves independently with the 
rigid environmental coupling. A first neighbour model in cubic 
Bravais crystals contains then only the parameters of 1 neigh- 
bour coupling (2 for sc < 100 > and BCC <111=' : f^^, 3 for FCC 

<110 > : fj^#f^,,f^) while is given by (1,36<): 


For SC: 


0 


[lOO] 


0 

0 


0 

3 

0 


0 

0 

3 


f^ = a ^ ^ 


J 


Stud 





(o) 


X 




(h) 

tr (12 ) 


2(fj^ + ^ 


(1.37a) 


For BCC; 


0 


llllj 


c; e C 

i' a 3 

0 r< a 


. £^=0 + 28/ = a - e 


and 


(o) _ 8 


3(fj_ + 2f^) 


(1.37b) 


For PCC: 


0 


[110] 


a 3 0 
0 a 0 
0 0 Y 


and 


(o) 


; £^ = a+0 , = a - 0 , = Y 


4 (f , + f , + f . ) 


(1.37c) 


S't. 

Figure (l,3e) illustrates the 1 neighbour model for the PCC 
lattice , 

Now wo discuss the translational-rotational invariance 

(h) 


I 

h 


(h) (h) 

0 and Z X * ~ 

B “ 


0 


(1.38) 


Validity in substitutionally disordered lattices. For simplicity^ I 

I 

we consider a sxxbstitutional defect/ as shown in Figxire (1,4a)/ j 

{ 

which disturbs the translational symmetry of the lattice. Now I 

y (a) > ! 

£ 0 0 0, Moreover/ the transversal springs between two atcstiS/j 

h j 

In general/ do not really represent an interaction between those j 
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(a) 

subatitutlonal 
represent ation . 


(Cj 

" spring ii ander tension 

tr , <&) Cutting z s^jring requires forces 
«V'ln order to ;!>ii.in'caia equilibrium 
Tb) far a transversal displacement , ip 
oi atom ^ tb® net force on atom 0 is 
tr'f jf tbifc caT« bo represented by a trans- 
voreal tpring^ wv~vf/ya « oy^. 



Cd) («) (f) 

Violation of rotational invariance for a aubatitutionally infect 
siociel (8C atroctura with lattice distance a), (a) Displacement 
pattern for rotation by an angla ''f about atom J. The forces on 
1 by the aprinji {0»1) — ftY * ft^* cancel* 

The same holds for th« forces from 3 and 3. h defect wltii a 

substituted at 0, The forces on one due to 

, i) Force pattern arou^ defact due to 
rotatiooi the resulting torque is 


different speihg fti 
0 and 2 are not balanced 


F(G.1*4 



two atoms alone but rather many body interactions. This is mor^ 
easily seen for a longitudinal spring f under tension o v;hich 
produces a ^transversal spring =cr/a (Figure 1.4b). The tension 
is carried by the environment. The consequence is that the modi- 
fied transverse springs violate rotational invariance in the 
lattice as shown in Figure (1.4c) for a SC structxare, i.e. total 
torque is non-zero. This simple example demonstrates that the 
springs - transverse and longitudinal, have to be properly treated 
if one wishes to study the substitutionally disordered lattices 
with force -const ant disorder (Chapters II and III). 

So far we have discussed the highest symmetry directions 
in real space which fix the eigenvectors and the degeneracy of 
tho eigenvalues for the coupling. The same considerations apply 
to ^(k.) in K-space if k is parallel to a high symmetry direction. 
Now we make use of these symmetries in defining Green's functions 
and dispersion curves. 


l.S Green's Fionction and Dispersion Curves 

The translational symmetry shows further that if e(ii) 
mn 

eigen vector of so is 


T(h) f e (n ) > 


e(n+h) (1.39| 

with the same eigenvalue, which means that by shifting a displa- | 
cement pattern its property of being an eigenvector of iZf(h) j 

I 

remains unchanged. This allows us to find immediately harmonic | 
solutions of the equations of motion by employing displacements ^ 
which are eigenvectors of T(h), namely, { 


a(n) *= A e 


ik.R 


,n 


ik.R^ 

e "" e(n) 




•fAfn) } 


e(n + h) 



with an eigenvalue exp Cik.R^) and with a yet unspecified "polar- 
ized" vector A. Thus 

/ \ 

-ik.R 


„ (h) ^ik.R*'-^ 

Ae - 1 ^ e‘ " Ae 

h 

•1 

= ) 2 f(k) A e 


.(m) 


(1.41) 


where 


- „ (h) “ik.R 

)2(k) = Z e * 

h 


(h) 


(1.41a) 


^(- k) = ^*C]s) 


has the properties 

^(k) 

Therefore the eigenvectors of j2! are determined by the 
eigenvectors e(k/cr)y a ~ 1/ 2/ 3 of the 3x3 symmetrical matrix 
0(k) [m ^^(k) = 0(k)]. 

j2f(k) e(ky c ) = ^(k' ^ ^ ) " 


with 

(e(k/ o ) y e(k ' , o')) 
Finally/ 


^ , e (k - k' ) 

o o 


(1*424 


^lka> = ji(ka)|ka> / <a. ]ka > = ~ 

/Vm 


ik.R 


m 


(l,42i 




All quantities are periodic in the reciprocal lattice 
i.e. k can be restricted to the Brillouin Zone/ k in V^. Thus 


(h) 


0 


(h) 


/ 


dk 


ik.R 


V 

Vb 


j2(k,) e 


(1.43| 


(in) ^ 

Since G ■” (w) and G(k/ w) also have same symmetries like 0 ancj 


^(k) respectively/ we can define/ 


i 



(1.44a) 


(h) 

(h) dk ijS'B 

G (w) = / “ . G(k, w) e 


o-. (k c ) e . (k 0 ) 

2 (1.44b) 

^(k ) - M(w + i n ) 


Here j2(k "^ ) are the eigenvalues of ^(k) 


^(kfT) = Mn^(ka) = (e(k .-), ^(k) e(ka )) 

which represent twice the potential energy per atom for a wave 
m Si 

S = e(k cr) exp(ik.R ). For cvibic crystals we use the results 

discussed in Section 1.4,1 where the eigenvectors of j3(R ) were 
h 

derived for R parallel to the main symmetry directions. The 

h 

eigenvectors of j2f(R ) illustrated in Figures (1.3a-d) can be 
viewed as polarizations for the corresponding three k-directions. j 
Dynamical stability requires i2(k) > 0 except for translations, 

* fm 

k = 0, = 0. The eigenvectors | k <? > form a ccmiplete basis for i 

I 

displacements with (S, S) < (x> ; the eigenvalues themselves cannot ; 
be normalized, but their displacements r>-'main atleast bounded at ; 
infinity. This does not hold for rotations. Consequently, rota- ; 
tions are not even admitted as proper displacement fields ; 

(Leibfried and Breuor, 1978). I 

To avoid the difficulties connected with the normaliza- 
tion of eigenvectors and the invariance inconsistency, we reverse | 

* I 

the proced'ure to get the dispersion curves in disordered lattices | 

I 

(Chapter III). If G(w) is known, we can get immediately G(k, w). | 
The imaginary part of G(k, w) directly gives information about j 

the eigenvalues. In a periodic system it is a delta fxinction at i 

[ 

eigenvalues, while a disordered lattice broadens the peaks making j 


V 


B 


with 


G. .(k, w) = I 

J a 




the phonon life time finite. For periodic systems we still use 
force constants model for evaluating the dispersion curves which 
is discussed below. 


1-5.1 Born- von Karman Models in Ciibic Lattices 

Because a theory of the potential and of the force para- 
meters is rather complicated and it is difficult to obtain results 
from first principles ^ one often uses models with fitting para- 
meters to represent the dispersion curves. The number of parametei 
should be small. In metals, as a rule, the electrons screen the 
interaction between ions such that one expects interaction only 
between near neighbours. Consequently one can use near neighboxar j 
force constants as fitting parameters. This is called Born-von I 
Karman model. If fi(k) is smooth, only near neighboxars contribute j 
whereas structxire indicates contributions from more distant neigh- 


boxirs. Here we restrict ourselves to short range interactions 
and for demonstration we choose the FCC lattice with 1 neighbox:^ 
interaction only, Pigxire 1.5, for which, | 


h 


(h) sin^(k.R~) 


4y (sin‘^ a + sin"^ — a) 


^ k. + k-5 ^ k. - k<, „ k- + k- 

+ 4 ct (sin^ T a + sin*^ — — r a + sin'^ 


+ sin 


2 


k. 


a) 


4a = 2(f^ + f^,) , 


4 Y = 4f , 




k k. -^Ic 

sin^ a), 4P = 2(fj,-f^,), 


1 . 



Comparison is made with experimental points for Ni. How much 
the deviations affect the study of disordered lattices will be 
treated in detail in Chapter III. 

This model is ratner abstract/ in the sense that one 
does not see the cause of the interaction. In particular/ for 
transversal springs/ tha origin is not clear and one cannot remove 
them without violuting rotational invariance. Consequently/ the 
use of such model for disordered lattice calculations is rather 
restricted. One can howover/ explore other couplings which are 
based on more physical models and which, then must obey rotational 
invariance. But we will not enter into a discussion of sophis- 
ticated models which are out of scope of our work. 

We have so far discussed Green's functions and their 
relations with dispersion curves. But there are also other exper- 
imental quantities viz. density of states, neutron scattering 
cross section, etc. which are related to Green's functions and 
are the check for the study of ordered and disordered lattices. 

So we discuss them for the sake of complateness. 


1.6,1 Density of states 

The vibrational (phonon) density of states p(w) is 
related to the "mass weighted" displacement-displacement Green's 
function G(w) (Eq. 1.21), 


p(w) Tr (Im M G(w 4 - i S ) ) (1,46) 

irdn 5 1 

-40 j 

{ 

where M is the mass matrix, n is the number of atoms per mit I 

I 

cell, d is the dlmensicaiality of the system. For the periodic. ‘ 



lattices, we can make use of Eqs, (1.44) to rewrite the density 
of states in terms of ^ (2s # 

p(w) = I ", (kcf)) (1.46a) 

jj. a 

This expression holds good if one gets the normal 
modes easily (Maradudin et al, 1971). But in disordered lattices, 
the transformations like (1.4 2) seem to be quite difficult because 
of lack of translational invariance. So we avoid ’’k-space" idea 
in our work wherever it is possible. Instead, we choose the real 
space technique - 'Recursion' method for evaluating G(w) and 
developing our formalism for the study of disordered alloys 
(Chapter II)- The recursion method is explained in Appendix A. 


1.6.2 Inelastic Neutron Scattering 

Inelastic neutron scattering cross-section consists of 
a coherent and an incohejrent part (Elliott and Leath, 1975; 
Sjolander, 1964). Tho coherent part is related to a weighted 
average over the displacement correlation function, specif i«hlly 


d^ cr ^ ^ h ^ 

d fidw ~ K 


r lx AA ‘ 

I q„. Ini <G^, (q, w)> n(w) 

aa» 

AA' 


(1.47) 


where 


B^(l) = ^>^(1) OXp{-'|<[2 . U^(1)J^>} 

“ (exp (-hw/KgT) - l y 



Here A defines the atomic species at site . The 


incoherent scattering 


^2 jji 

d fidw 


h 

IT 


K 


I 

aA 


Im. <G^ (0, 0; w)>n(w) 


(1.48) 


is related to a weighted average of the single site Green's 
fxmqtion and hence to the density of states. 

These weighted Green's functions are related to unweighted 
Green's fionctions. For single site approximation the ex;pressions 
developed by Elliott and Leath (1975) are frequently used. But 
we follow a entirely different approach. After obtaining confi- 
gurationally averaged Green's functions (discussed in Chapter II )y 
we make use of 'Recursion' method (Appendix A) for getting 
weighted Green ' s functions . Here we carry out the similar recur- 
sion procedure as in the ordered lattice with new self-energies 
except at those sites where we wish to evaluate the weighted 

i 

Green's function. j 

i 

However, th':^ coherent cross-section allows a direct ' j 
measurement of <G(k, w)"^ and gives the most detailed experimental ! 
information about phonons in mixed crystals. The imaginary part | 
of <G(k:/w)>, the spectral function, is proportional to line- | 

shapes and their widths are finite. These are the quantities | 

I 

which are frequently measured and hence we discuss them in | 

I 

Chapter III . [ 



Chapter II 


The Coherent ?otential Approximations 

2.1 Introduction 

In this chapter we shall review briefly the existing 
methods for the solution of Eq. (1.24) alongwith our method. 
Though attractive at firsts some approaches become intractable 
for the phonon study, while others involving simplifying assump- 
tions either become far removed from reality or lob^e features 
essential for physical relevance. In the study of phonons in 
disordered solids all the disorders - diagonal, off-diagonal, and 
environmental - are coupled. Our aim, theiref ore , is to examine 
these methods with their shortcomings sind discuss the basic ideas 
behind the different approaches, so as to expose avenues of 
generalisation including all effects which yield physically 
valid resiilts in the real alloys. 

2.2 Virtual Crystal Approximation (VGA) 

I 

f 

In the virtual crystal limit when the scattering by the J 

disorder is very weak i.e. the mass changes and force constant 

_ ,, i 

changes are small, a perturbation expansion in powers of V as ; 
defined in (1.24) is valid. The expansion is 

f 

<G(1,1«)^ = G°(l,l') +I G^(l,m) <v(m,m' ) > G°(m‘,l* ) 

iraTTi* 

+ I G°(l,m)< y(m>m*) G°(m',n) y(n,n') >G°Cn‘,l*J 
Iran' 

nn' (2. 1 





In the VCA^ it is ass\araed that all tho*V(mym') are very small 
so that the higher order products of V's need only be treated 
approximately. Thus all such average products of V's are decoup- 
led (a random phase approximation) into products of average 

<y(l,l') V(2,2') ... V(n,n')> 

y<V(l,l' ) ><v(2,2' ) > ... <V(n,n')> (2.2) 

This gives 

<G(w) > = G°{w) fl -<y >G°(w)]'’^ (2.3a) 

= (2.3b) 

so that the Green ' s f lonction is approximated by one for a perfect 
crystal with the atomic masses equal to the mean mass of the 
disordered crystal and with force constants equal to the mean 
force constant of the disordered crystal. Thus^ the frequency 
spectrum of the disordered crystal is rigidly shifted to that of 
a virtual crystal without broadening the states. The phonon 
dispersion curves in the VGA limit simply shift linearly with 
concentration. Here phonon K-modes have on infinite life time, 

2.3 The Average T-matrix Approximation (ATA) 

in 

Rearrange equation (2.1)ysuch a way that all terros invol- 
ving repeated scattering by the same impurity are collected 
together/ 

<G(1/1»)> = G°(l/1') + I G°(l/ra) < t(m/m‘ ) > G°(m' /1‘ ) 

m 

+ 1 ' G^d/itv) <t(m/m.') G°(m'/i^) t(n/n')> G°(n'/1‘) 
m/n ” 


(2.4) 


where t(m/m') 


. 1 ^ . . 

= t IS given by 

t = V (1 ~ V)'^ (2.5) 

In matrix form 

< G > = 0*^+0° <t > G° + <t G^'' t > G° + ... (2.6) 

where 

G*^ (n^m) = G°(n,m) ■” S nm‘ (2.7) 

The average t-matrix approximation (KTh) is to decouple the 
average of the product of t ' s according to 


t(2) ... t(n)> = <t(l)><t(2)> ... <t(n) > (2.8) 


which is somewhat better approximation than the VGA (2.2), 

o * 

because adjacent sites are restricted by G from coinciding 
(Elliott and Taylor, 1967). Then, 

<G > = 2° + G° J< t > (I - G° < t > )'^] G° (2.9) 

or ^ as defined in Dyson's equation 

"^2^ ~ 2*^ 4- G^ <G> is given by 


£ == Z;'! + <t> G^'J (2.10) I 

j 

For mass disoruer V is diagonal, in the site ire presentation and j 

I 

I 

V.„,(l,l') = AM„w^ « 6(1,1') (2.11)1 


and so is t. Then 


< t > 
aa 


C AMjj W 


|i — A w 


( 2.12 


and 



(•2.13) 


Z (w) 
cia' ' 


C A Mq W 


[l - (1 - c) AM„ G°^(o) 


This single site self energy is site diagonal and is 
eqxxi valent to an effective mass (complex and frequency dependent) 
on the disordered site. A weak point of this approximation is 
its general failure to locate the band edges . To shift the band 
edges from those of the unperturbed lattice one must treat the 
scattering in a self consistent field. 

2.4 Single Site Mass Coherent Potential Approximation 
(ICPA or MCPA) 

It is a statistically averaged theory which specifies 
an average Green's function for the equation (1.24) (Taylor/ 

1967; Soven/ 1967). This average -medium Green's function £ is 
fixed by the requirements (i) that the effective -medium quasi- 
particles scatter from each atomic site the minimum amount/ that 
iS/ the single -site effective medium transition matrix/ when 
averaged over all possible alloy configuration is zero (it is 
not practical to require that the multisite T-matrix vanishes ) , 
and (ii) that its self-energy assumes the mathematically 

simple form I 


I ln> cT_(w) <nl ^ I 
n n 


(2.14 


thereby giving a single frequency dependent life-time and level j 
shift to all the quasinormal modes of the statistically averaged ! 

' i 

alloy. Requirement (ii) gives the simple result for the effecti-j 
mediuim Green ‘ s function i 


£(w'^) s= G°(w^ - m~^ cf(w)) 


( 2 . 1 ! 



o i2> 

where G°(w'^) is the perfect lattice Green’s function for a lat- 
tice with all masses equal to rag. In the mCPA-^ the choice of 
reference lattice is irrelevant, because the self-consistent MCPA 
Green's function is invariant under a change of reference lattice 
(Elliott et al, 1974 ). Rewrite the equation (1.24) 

G = G° + V G (2.ie> 


where we have 

V = (g - ZJ - Cm - (2.16. 

For the mass-disorder only, V is always diagonal and is non-zero 
only at impurity atom sites: 


V 


with 



(2.17a) 


-n 


2 

I n> (rOg - itij^) w'^ <n 1 


(2.17b) 


The effective medium Green's function is defined as the confi- 
guration average of the alloy Green's function G, 

£ = « G » (2.18) 

where double braclets denote an ensemble average over alloy^-ocaa- 
figurations. In the M'PA the statistically averaged Green's 
function is assximed to satisfy 


£ = + G^ L. 2. (2.19) 

which results in the simple expression Eq. (2.15) where £ is 
the (as yet lanknown) self -energy. Eliminating G° from Eqs. 

(2. Id) and (2.19), we obtain. 



S. = 2+£(Y- £)£ 


a + ^ - c ) G 

•* ^ n n — 


( 2 . 20 ) 


The single -site effective modi uni transition matrix is 


(v - a ) 
n n' 


n 


fl - g(v^ 


( - - ^ V ^ a ^ . 


( 2 . 21 ) 


In the 1.CPA (mcpa), the configuration-average of must vanish ^ 
thus we have 


« T » 

n 


0 


n > 


< n I ( ■ L3: ... . , r **■ ■ £ ,, i ! ■ 


(lUg - m^)w^ - aj 


) 


1 + rr ^nig^n* ^ _ <n|g|n> [(rn^ - ra^)w^ - a] 

(2.22a) 

which results in the equation for cr(w)^ 


a(w) 


c(rag - m^)w'‘ 


1 -<nl G°(w^ - c(w)) [n > pm^ - 


(2.22b) i 

When combined with the known expression for G°, Eq, (2.22b) can | 

/ V I 

be solved for the self-consistent self -energy a(w). The self- j 

energy is diagonal and branch independent. The MCPA phonon | 

j 

density of states ; takes the form. ; 


p(w) 


Im (Tr « M G>>) 

N IT — 


(~) lm(rn^(l - 5 ) < n | G°(w^ (1 -5 )) |n>) 


o , 2 


with 


2 *" / \ 
Itlg W <7(w) 


(2.23a) 


(2.23b) 


^ (w) 



2.4,1 Failures of mcpa 


The MCPA provides a reasonably good description of an 
alloy over a wide range of parameters, it works well when mass 
difference of two components of* a binary alloy is small and 
becomes inadequate when . y .,.. ±s large. (Here Z is the 

F I 

number of nearest neighbours). 

Since MCPA is a single site approximation it cannot take 
correctly into accoxont the effects due to two or more sites. 

Firstly/ it cannot correctly account for the randomness in the 
force constant matrices involving two or more sites. In binary 
alloys the shape of the vibrational density of states function 
as well as the spectrum width depend on relative values of various i 

I 

interactions between neighbovuring atoms. These nuihbeirs «an be i 

very differen+- for the two constituents of an alloy even though ; 

their atomic force constants are similar. Secondly, short range j 

[ 

order due to chemical clustering cannot be taken into account , I 

Finally, the effect of statistical clustering which is parti- I 

i 

cularly important in the impurity bands is totally absent. The 

f 

local modes show some fine structures (Payton and Visscher, 1967| | 

Dean, 1972; Alben et al, 1975) which are lost in MCPA. The MCPA | 

fails to give correct bandwidths oven in the absence of off- | 

diagonal disorder. In fact, MCPA (ICPA) was not designed to j 

I 

reproduce the spectrum, but was constructed to properly mimic in 
long wavelength properties of the lattice vibrations in the alloy. 
The experiments on the vibrational properties of Ni-Pd, 

I 

Ni-Cr, Ni-Pt alloys indicate that local environment strongly | 

f 

affects the density of states, dispersion curves, and width ! 

functions (Kamitakaha and Brockhouse, 1974; Bosi et al, 1978; 1980; 



Tsunodct/ 1979), In fact WCPA ignores all environmental effects. 

The failure of ICP/^ (MCPA) for large y (defined in (1.24) can be 
understood from the fact that for large V the excitation mean 
free path is small and in order to obtain a good estimate of the 
configurationally averaged Green's function wo must consider 
scattering due to pairs ^ triplets etc. in a self -consistent manner. 

2.5 Generalisation of CPh 

Prom time to time suggestions have been put forward about 
ways of generalising the CP7i to taKe into account the effect of 
clusters, off-diagonal disorder and SRO,for the multisite scatt- 
ering, neglected in ordinary CPA (Mookerjee, 1979; Srivastava and 
Joshi, 1973). In principle, these generalisations should prodijce 
accixrate alloy theories, provided a sufficient number of sites 
are included in the determination of a self-consistent self- 
energy £ (k, w). ! 

In the multisite scattering approach the system is parti- , 
tioned into clusters, each containing n— sites (Tsukada, 1972) ; 

I 

and clusters are labelled as ^ 2 ' '** *** indivi— [ 

dual sites are labelled as R^/ R 2 / ••• Rj •• • Parallel to | 

single site transition matrix (2.21) | 

I 

•Lj = CVj - (2.24) j 

where Gj is an (n x n) matrix. The CCPA equation then reads j 

^ V - n (2.25) I 

<T^j > = U I 

The self-consistency condition (2.25) may be written in | 

terms of the matrix elements of the configurationally averaged 
Green's fxinction/ if we observe that the Green's function for the | 

f 

[ 
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system with the averaged medium everywhere except the cluster 
(or cell) can be written ass 



Then^ 


<G> +<G>Tj 


<G > 


(2.26) 


^ 2 j ^ ^ 2 ^ (2.26a) 

which is similar to MCPA condition. 

This equation when applied to simple model of random 
binary alloy on a linear chain, with diagonal-disorder only, yields 
satisfactory results . The band edges are produced more accurately 
than in MCPA. Ducastelle has shown that for CCPA with cluster- 
diagonal self -energy, the approximate Green's function retains 
the Harglotz property of the exact Green's function which is a very 
favovirable point for the generalisation in this form. But when 
we apply CCPA equations to realistic situations the calculations 
become infeasible (Tsukada, 1972), For this reason several other 
simplifications have be.en suggested e.g. central site approxi- 

I 

^ I 

mation, boundary site approximation etc. (Kumar and Joshi, 1975), | 

So far we have discussed the situations where only ! 

i 

diagonal-disorder was present. Several attempts have been made | 

' I 

to include randcxnness in the hopping integrals especially for i 

I 

electronic studies (Kumar et al, 1975). Here it is assurtKSd that i 
the diagonal term of the Hamiltonian can have either of two values i 

* f 

itu or and hopping integrals may have values 0^, 0^ or 0^^ j 

o A I 

depending on the configuration of the nearest neighboxir pairs. | 

f 

! 



2,6 Cluster EmhCiddinci 


In addition to various approaches for cluster general- 
isations discussed in the previous section ^ people have tried to 
embed largo size clusters in some non -sc If -con sis tent effective 
medium e.g. immersion of a large size cluster in MCPA or VGA 
medium (Diehl and leath, 1979a; Myles and Dow, 1979; Choudhry, 
1981). Here once the effective medium (MCPA or VCh) is obtained, 
we get the effective Green's function ^(w). Embedded in the 
effective medium is a cluster of N atoms giving the Green's 
function for the alloy with this specific configuration of atoms 
in the cluster G(w), We then define the scattering potential 
V*(w) which vanishes outside the cluster. The alloy Green's 
fmetion G can be written as 

G = 2. 2. X 2. (2.27) 

which need only be solved for atcans within the cluster (G = £ : 

I 

outside ) . This gives i 

( 

G = (^ - 2 V ) "^ 2 (2.28)! 


By choosing cluster of larger sizes, details of snvirc^i- 
ment implanted can bo gradually increased. For each possible | 

.configuration over a cluster, one can determine the density of ! 

f 

states at the central site and also the configurationally averaged | 
Green's function, particular features in the density of states | 

‘ ' f 

curve can be attributed to specific clusters . j 

But this kind of cluster embedding often gives rise to 
two kinds of difficulties: (-1) the band width remains incorrect 
and (2) because the density of states per mode still integrates | 



to unity, this has the effect of distorting the structure speci- 
ally near the band '.jdges to preserve normalization. Thus spurious 
structure may be obtained where there should be none. 


From the 'discussion of the previous few sections, it is 
obvious that none of the proposed cluster generalisations are 
satisfactory. We, therefore, turn to on altogether different 
viewpoint of configuration-averaging, namely the "augmented space" 
(Mookerjee, 1973; Mookerjee and Bishop, 1974; Haydock and 
Mookerjee, 1974; Mookerjee, 1975a; 1975b; Kaplan and Gray, 1976-77 
Kaplan et al, 1980; Kumar et al, 1982). Basing ourselves on this 
formalism, we shall develop an ijinambiguous , tractable general- 
isation of lCP7i (MCPA) which retains the Herglotz property at all 
frequtencies for all ranges of disorder. Off-diagonal disorder 
and short range order shall also be incorporated. 

Recently in a series of papers (Kaplan and Gray, 1976; 
1977| Diehl and Leath, 1979a; 1979b) the authors have examined 
"augmented space" and have generated self-consistent resxilts, 
with off-diagonal disorder on linear chains in particular (Kaplan 
et al, 1980), Parallel to their- algebraic equations is the 
graphic al method (Mookerjee and Bishop, 1974) that has been 
adapted in our work. The graphical method has the added advan- 
tage that it is ideally suited in conjunction with the "Recxirsion" 
technique (Appendix A) of Haydock et al. (1972, 1975) for the 
application to realistic models on three-dimensional lattices. 

The equivalent algebraic equations soon become tediously lanwiel- 
df, as soon as the cluster size increases, particularly on three- 
dimensional lattices. 


/ 


I' 

1 


r 



2.7 ' Augmented Space ‘ Formalism 


The Green's function. (1.21) or (1.24) for the disordered 
binary alloy's Hamiltonian (1-2) is averaged in the earlier 
Sections (Sections 2,1 to 2.6) where basic principle is same: 

The dynamical behaviour of the system is described by a Hamiltonian 
suitably set'up^ whereas the statistical behaviour is^ as it were^ 
imposed from the outside by allowing the potentials involved in 
the Hamiltonian to vary randomly according to some prescription. 

The Hamiltonian/ in itself/ does not describe the full behaviour 
of the system completely and has to be augmented with the distri- 
bution of the set of random potentials (the conf igiiration ) • In 
the 'augmented space' formalism, the excitations in a random array 
of disordered potentials are considered to be in periodic poten- 
tials aiong with a "field" which describes the disorder (Mookerjee, 
1979). An analogy is the case of electrons in the presence of an 
electromagnetic field. The Hamiltonian, expanded no include the 
disorder field, completely describes the disordered system. It 
is not necessary to put in the statistical description by hand- 

Sinco the statistical description is incorporated within the ! 

I 

Hamiltonian, the configuration averaging is not a further extra | 
process as in the earlier mentioned approaches. ' 

First we shall give a version of the ‘augmented space* ' 

t 

formalism (ASF ) of Mooker jee for the disordered vibrational system | 

case that we are interested in. We start with the assumption I 

! 

that the lattice sites are randomly occupied by the atoms of type j 

i 

A and B in a binary system. The Hamiltonian involves a random | 
occupation variable defined to take the value 0 or 1 I 







corresponding to the occupation of site i by A or 3 type of atom. 
A set of values taken up by the random occupation variables 'f 
is called a configuration. Suppose there is no correlation 
between these random variables so that we can define their distri- 
bution f\inction in the form 


p (^N^>) = n 


(2.29) 


where is the probability density of the individual vari- 

ables. This assumption ignores short range order due to chemical 
cluster effect. The more general case of dependent variables 


have been considered by Kaplan and Gray (1981). 

Haydock (1972) noted the fact that the probability density 
Pi(N^) satisfies the propertiejs 


Pi(Np 


p^(N^) dN^ 


These properties arc specific to the imaginary part of a Herglotz 
fionction e.g. the resolvent of a self-adjoint operator in a 

certain Hilbert space. Taking clue from this Hooker jee (1973) 
suggested that a Hilbert space '-f and a basis {| f^ >} may be 
introduced such that P^(N^) corresponds to the imaginary part of 
the resolvent of a suitably chosen operator in If lfQ> 

is a specifically chosen member of the orthonormal basis ^ I f ^ 

in then I 

I 

I (2.30) I 

where Z -»■ N_i + i 0^ ! 

i 

For a given p(2), one can construct a suitable operator j 
M provided p(Z) can be expanded in a convergent continued frac- ! 



41 



2 - a3 “ 

m 

m 

m 

The operator M v/ill then be a tri diagonal matrix with ^ 2 ' 
a^ ... etc. along the diagonal and 102 * ... ... etc- along 

the off-diagonal positions. 



in some basis ^ convergent continued fraction 

expansion / we have to restrict oxarselves to probability functions 
with moments finite to all orders. This conforms to physical 
requirement. Most of the physically valid probability distri- 
butions are of this kind/ the exception only being the Lorent- 
zian. This is the inverse of the procedure involved for deter- 
mining matrix elements of Green’s operator from the recursion 
method ( Haydock / 1972). 

For a binary alloy A B. with diagonal disorder only 
{ mj^ } forms a set of random variables described by 

“i “ ny + rog (1 - N^) 

where ^ } are a set of random occupation variables for the 

solute A in the solvent B. is a kind of indicator function 
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which carries all the conf igxiration dependence information. In 
the absence of short-range order the various are independent 

random variables v^ith probability densities. 


p(N^) 


c ' (N, - 1) + (1 - c) 5 (N, ) 


(2.31) 


i - r C 1 - C 

^ N + i ^ _ 1 i-7 + i ^ 


Im 

4- 

n-vO^ 


— + 1 ^ N ‘ = N . + i n 

N* - 1 N' ^ 1- j- • 


= - T Im { 


ri 0 


+ N* - c 


c(l - c) 

N* -TI”- c) 


Im ( 
n-^O' 


(2.31) 


N' 


1 N' - a. 


Where 


c, 3-2 = (1 - c), = (1 - c)c. 


A representation 


of the matrix M is thus a 2 x 2 matrix 

0 

c ’^cd - c) 

c(l - c) 1 - c 


/ 


(2.32) 


with eigenvalues 0 and 1- Its eigenstates and 1^2^ 

describe occupancies corresponding to these eigenvalues. 

We now introduce the "disorder field" . To each variable 

N. we have associated a vector space . and a representation of 
X 

a self-adjoint operator corresponding to the probability 
density of that variable in some basis {] f^>}. If we had 
taken the set of eigenvectors {] h^ >} of as our basis then 
each 1 > in would correspond to one particular value ta3cen 



by the variable In general any element in is a linear 

combination of the {] «' 'f j_ with its basis and operator 

thus contains all possible informations about the configura- 
tions of the variable 

The product space ^ y, can now be constructed 

i 

which contains all possible states of the set {Nj^} i.e. all 
configurations of the disordered system. The basis in the 
product space is 

l£> = {lfi>} ® ig) . . . (2.33) 

(Si.) _1 

and if 2 (z) = (Zl-Mj_) 

then 

p({N }) = Im<f lG(z)!f > 

j« tr vj w 

where 

|fo> = ® 

and 

(Mt) - ^ (M3) ^ 

G(z) =2 ^ (Z) @ 2 (Z) ® 2 (Z) ® . . . 

We call 'i' the disorder field, it is clear from th« way it is 
constructed that it contains complete information about the 
statistical behaviour of the system. The information about 
probability of configurations is described by the operators 
defined on the field. 

We now describe the configuration averaging procedure. 
Let us first consider a function f (s) of one variable only. Its 


average is 


C30 


f = / f(s) p(s) (is 

- CO 


.ao 


= / £ (s ) { - 


CD 


Im 

Z s + iO 


+ 


(Z) } cIS 

-“oo 


Assuming now that £(Z) is a fxmction of a complex variable Z 

and has no singularities on the real axis, in the neighboxarhood 

(m) 

of the branch cut of the function g^— (Z), 


1 


2 TT i 


f(z) gQ§^(z) dz 


(2.35) 


The contour over which the integration is carried out is taken 
into account the branch cut of g^— ^(Z) along the real axis and 
not including any singularities of f(Z). Since M is a Hermitian 
operator on 'f we can write , 


00 


^ £(Z) dZ <f I f (Z - h)~^ dp(h) |f > 


2 IT i 


- CD 


where we have written ^^—^(Z) = / (Z-h) ■^d_ 2 (h)/j 9 (h) being 

- 00 

the spectral projection operator of M. Thus, 

— ^ 

£ = ^ lfQ> 

-CD 


CD 


-1 


'^o^ 


(2.36) 


where f (M) is the same function of M as f is of s. The average 
has been expressed as a representation of a suitably constructed 
operator . 

The generalisation to a function of several variables 
is straightforward with caie associated with each N^. In case 
f(N) is itself a matrix element of the kind say <r^| F(s) 1 r ^ > / 



F £ i#tato vectors jr^j f^> = & 

the augmented product space Q ~ Jt X '¥ 

F = < f o I / ^ G^~^(h + iO"^) } 

- OD 


! f > are defined on 
' o 

and the average becomes 


F(h) dh I 5,^. > 


(2.37) 


an matrix element in the expandea space Q . 

Now we are in position to use 'augmented space' in the 
vibrational stud^ii' of siobstitutionally disordered binary alloys 
which is discussed in the following section. 


2.8 Lattice Green's Function in 'Augmented Space' 

We wish to calculate the conf ig\nration averaged one- 
phonon Green ' s function 

<G^j(w)> = JJ ... / G^j w) n p(N^) dN_^ (2.38) 

whore the displacement -displacement Green's function, defined 
in (1.20), for a binary harmonic alloy with nearest ' 

neighbour forces is | 

i 

Gij({Ni}, w) = <r^l 1 rj > (2,39a) [ 

I 

with 

+ Zl .(N. , N ) T ? (2.39b) | 

I 

Here 1^^^ basis for the Hilbert space and the Integra- | 

tion over p(N. ) dN. constitutes averaging overall possible confi- 

Ju JL 

gurations of the disordered system. = I i > < i I and 




I 


(m(N^) w" 


-I 


l^i 




:il(Ni, 


N, 



T^j = !■!*■< j the- projection and transfer operators in the 

Hilbert space- on the tight -binding basis. W„ have incorporated 

in (2.39b) the Siam rule 1 0. .(N., N.) = 0 essential to climi- 

j i' 3 

natc the w = 0 mode corresponding to boaily translations of the 
system as a whole (conservation of momentum). This siom rule 
was earlier cumbersome to incorporate in the traditional CPA 
calculations. The consequence is uhe environmental disorder 
which couplos diagonal and off-oiagonal disorder in the Hamil- 
toaiian g (second term on the right hand side of eq. (2.39b)). 

As a result/ the aiagonal element ^ / and hence G ^. , depends 

on the occupation of neighbouring sites. { mXN^^) } and 

{^. j(N. , N.) } are the sot of random variables described by 

Xj X J 

m(Njj_) = I + m^d - N^) I 

^ij(Ni/ N^) = ^ Nj 4- - Nj_)(l - Nj) 

+ {N.(l - N.) + N.(1 (2.40) 

where arc a sot of random 'occupation' variables for the 

solute h in the solvent B. In the absence of SRO the various 
Nj^'s are independent random, variables with probability density 

p(N^) « c 5(N^ - 1) + (1 - c) 6(N^) (2.41) 


using 
6 “ 


(2.40), wo may rewrite (2.39b) in more useful form 

K„ + I mw^ I N. P. “ I (I 4- 

i - 1 i i jL^^i. 

4-^^^^ N_.N^) P^ 4- I I 4- Nj) ) T^. 

r j 


(2.4Z) 



where m = J - ^^7 0^^'^ = “ 2^ 

ana K^ is the non rcindom operator for the pure solvent B, 


vBB 


•4B 


K 


B 




I . ^BB ^ I ^BB _ 

— ^oo' —1 . . ^ — JLJ 

X2 


(2.42a) 


The opjrator K is now in a form appropriate for the 
‘augmented space' trensforrr.ation. Corresponding to each 

• (i) . ^ 'f (-i-) 

there as a operator M j.n a two damensaonal Vector space 

A representation of in a basis {iv^>^ |v^>}is 


»'cTi~~~~'cy 


^c(l ~ c) 


1 - c 


; (Mookerjeo, 1973) 


where c is tho concentration of the solute A. The augmentoa 
space operators, now, are 


K 


I w" - ,I 4i)£i + I lij 


li^i 


ij 


^ I + mw^ I I -P. 0 (X) I 


(i) 


i 101 


i J 


(±) 


-ij 

(2.43) 


K j(ij) the above equation indicate identity 

- - . 

operators in all subspaces J except those superscripted. 

The configiirationally averaged Green's function (2.38) 


is, then given by 



4b 


<Gj_j(w)> = <fj_ t'“' i i rj > (2..44) 

with j > = r > . 

i • 

4S«t ^ 

or in torms of m 1, 

^ < r j_ I (ra w^ i - ^)**^ 1 r^. > (2,44a) 


The cxpond^i operators m, ^ anu k contain complete information 
about the system and are defined on the expanded Hilbert space 

r ■ . 

Q, = ^'\‘x 'i' I'^here 'i' = n '(f . . Description of the dynamical 

i 

behaviour of the system is in the subspace , while tnat of 
the statistical behaviour is in 'f , the disorder space. The 
configuration averaged Green's function is a particular repre- 
sentation in this expanued space n . The problem has thus been 
reduced to the determination of the matrix element of the 
operator g in the augmented space. Notice that upto this stage 

of configuration averaging, no approximations are involved. 

¥ 

Subsequently/ in- th- determination of tne matrix element in the 
a\jgmented spiice herglotzicity retaining approximations may now 
be introduced systematically. In the next section we shall 
describe the Graphical mothod which will enable us to do just 
this . 


2.9 The Graphical Method 

The graphical method was introduced by Haydock (1972) 
for calculation of the Green’s function. It will bo used to 
generate the cluster CPA formalism for a general alloy system 
in Section 2,10. 



fly a- graph we shall mean a set of ‘vertices* connected 


by the 'links*, lat {li>^ be any complete linearly independent 
basis in which the Hamiltonian is described. The overlap matrix 


^ij 


<± ! j > is a unit matrix in case the basis is orthonormal. 


Matrix element of an opvjrator M are given by M. . = <i 1 M I j 

J 

A graph is associated with every basis {ji >> . To every element 
1 i > of the basis set { ! i > } is associatea a vertex v^ ana to 
each distinct pair of elem-w-nts (ji>/ ! j >) a link or bond 1- . 

which may be directional (i.e. l^j may not be equal to ljj_) 
otherwise. To any operator M there corresponds a graph consis- 
ting of vertices and those links for which M. . 0, If we want 

-i- J 

to invert an operator M, then the contribution of each vertex 
is defined as 


k(vp 


1/M. . and that of the link 1. . as 

X jl 


k(l. . ) = M. . 

ij 

A path Qj^ of length N is defined as a series of ordered 
vertices and links i.e. 




• ''n ^ 


Contribution of a path is defined as 


k(Qjg ) 


n 

vertices 


k(Vi) 


n 

links 


kdij) 


A self-avoiding path of length Qj^ “ ^ ’'^’i ^ij ' "^n^ 

a path such that v^^ = v^ with the only exception that v^ may be 
the sane as v^/ in which case it is a self -avoiding closed poly- 
gonal path* 





The renormalized contribution of a vertex v in a 

n 

path is dofiii.^d as 
00 

^ I ]c(Q«) (2,45) 

N=0 N 

where Q^l are all polygonal self-avoiding paths from v^^^ and back, 

on the graph from x-'j'hich the vortices v^ . . . have been 

-1 

removed. If H = G then the representation of G in some basis 
can be expressed (Heydock, 1972; Mookarjee, 1979) as follows; 



CD 

I 

N=0 



k(Qj,) 


(2.46) 


where are all self-avoiding paths from v^^ to Vj. 

Nov; we make use of this technique to calculate the 
Green's function defined in Eq, (2.44) as the technique as such 
presexves the Herglozicity (Haydock, 1972 ) . The essential 
approximation now involves ignoring or delinking the contribu- 
tion of possible self-avoiding paths. The relation of these 
polygonal paths to the CPA has been discussed by Bishop and 
Mookerjee (1974) and Mookerjee (1975b) has demonstrated their 
relation to the multiple scattering diagrams. The philosophy 
of approximation is as follows. 

We consider exactly all paths that involve only the 
spatial part of the augmented space i.e. >£ . Such paths are 
connected to the lattice structure and must not be approximated 
or delinked without destroying some of the feature of the under- 
lying lattice. In 2- and 3-dimensions the recursion method 
helps us to deal with this infinite set of paths. 



if we want a cluster-CPA with ir^/ ^ 2 * ••• within 
the cluster, then all polygonal paths that do not completely lie 
in the subspace and involve these vertices aiad none others 
must also be exactly acco-unted for. These are related to the 
multiple scattering within the cluster. 

Larger polygonal paths which involve sites within and 
without the cluster but which do not completely lie in the sub- 
space y-i must be delinkea. The type of delinking depends very 
much on the environment we wish to immerse our cluster in. 

Immersion in a VGA raedixom yields the travelling pseudo fermion 
approximation (generalization of the ATA) of Diehl and Death 
(1979a). If we wish to immerse in a ICPA mediumy we simply 
delink all larger polygonal paths into a Cayley tree (Bishop and 
Mookerjee, 1974). The result is not self-consistent and in 
particular the bandwidth is the same as that of the ICPA. If 
we wish to embed our cluster in a self -consistent cluster medium 
wo partition the lattice in Tsukada typ«, clusters, delink the 
environment itself into cluster paths and self -consistently 
calculate the »jffact. Wo shall describe this in detail below. 

2,10 CCPA Self-Consistent Madium 

Here we shall illustrate the generation of a self -consist 
modiixn in the CCPA for oxxc phonon problem. For electronic 
problem, it is already used extensively by Sri vast ava (1982). | 

To stort with we consider the 2CPA. The generalization of bigger 

clusters will be dealt with subsequently. i 

[ 

S' 

A grajii is constructed according to the delinking proce- | 
dxire described in the previous section which has got the structure 
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FIG.2‘1 (continued) 

(see next page for details of symbols) 



Dictionary of symbols used in Figure 2,1. 


^01 

^l(6) 

26o^5) 

2ol = > 

2oo f ^ ) 

26 o <®> 


0®® + nr. «(!) L 2 of(2) 
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fc(l-o) (a^) + c 


( 2 ) 


-01 


01 
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-0(1-0) + c 4o') 


cd-c) (i^« + c(l-o) 2 oo 0 

c(l-c) (2®® +.2^2) ^2). 

+2oi’ +0(1-0) 2^2)) 

■ flics ') (2^“ + (1-c) 2oo') 

''o(l-c> (Mo“ + (l-<^) > 

>'c(l-c) (mgw 2 I + g(l> + (l-c) ^ 0 ^) 
( 2 g® 4 2 (l-c) + ( 1 - 0)2 gU)) 


= 1 + i ^ So ) + =(” 1 + %’ + = aio ^> 

- (nigw2 I + 2^®) + ( (1-c) (nw2 I + o 2^0^ ) + 2oo’ ) 

= (mgw2 1 + 2^®) + (l-o)(nw2 l + 22^g^ + (1-c) 2 ^q^) 

= (mgw2 I + ^ (o(raw2 ^ ^ (i.c))2^o' + 2oV 1 







and HK. Once the rt^normalization is accounted for^ the octagon 


is effectively isolated. The renormalized vertices o (w) and 

— o 

links ^ founa as tollows. Let us divide the lattice as 

follows into two subspcices (1) an unrenormalized bond AE, which 

we shall call system 1. It has a Hamiltonian = Z (m^w^ 

(i) r i=A,E 

^oo i ^ /2f. . ( X ) T. . ( X ) v;here x is the direction vector 

ij=A/E ^ 

of bond AE. (2) A lattice which is tnc original lattice minus 

the link AB, ana in which all bonds and sites are renormalized 
• This has a Hamiltonian 

= I P. + I (w; X ) T. h(X ) (2.48) 

iA/E iyjA/E ^ 

and ( 3 ) a linking Hamiltonian 


K 


inti.. 


= I 


£l(w. 


X ) ) + 


3:^(x) 




T^(X) 


+ 


T ( 
“Ei^ 


X) (2.49) 


where x is the direction vector corresponding to the bond ij 


subscripted in transfer ox^::rator j 


-P 


If we now ui;a the notation X ^ to denote the inverse of 


the operator ^ in the subspace Y 

= il § il 

where = Lw^ I - i'" 

with El + I 2 “ = (2.50) 

Examination of the equation (2.50) immediately indicates 
that the effect of the rest of the lattice hanging on to the bond 
AE is to change the Hamiltonian to + a where the self- 

energy %. is given by 



a 

-o 


a 

— AA 


-EE 


I I^CAk) ziJ’CntA) 


k/iti eN 


(2.51) 


-AE ~ -EA -1^^) -- ^ I _ l^(Ak) Z_^(^) 


r Or. 


ke 
m e 


where denotes a nearest neighbour of m. and AE denotes the 

direction vector corresponding to bond AE. r is the inversion 
matrix with respect to bond /iS. 

(b) Let us now work on the renormalized, isolated octagon. 

7iS before let us divide the octagon into two siobsystems; a bond 
B (Of^ If), the rest of the octagon M. 
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ij Of, If 


K. . , K 
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ijyof ,1£ 
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I I iij + Kji) 

ISB jcM 


/vs before , we have 


2' 
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Eb 2 £b 


r 2 


[w" i - gB - 


BM ”BM-t -B 




» "[w^ I - K*'"J 

This immediately shows how the octagon renormalizes 
the bond B. We have 




H 


,’;bm _m -6m 
— oj — jk “^o 


(2.52) 


r- ! \ J. Z V tBM 

ll(AE) = 1^1 -jk ^1 


Equations (2.51) and (2,52) together provide a self- 
consistent set of equations for the calculations of the 
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self-energy ^ , The various Green functions involved' can be 
easily calculated on a three-dimensional lattice using the 
rec\arsion technique of Haydock et al. (see Appendix A). In the 
case when we are using frequency dependent effective Hamiltonians/ 
the recursion has to be done separately for each value of freq- 
uency. 

The cluster generalisation to arbitrary large clusters 
now follows exactly the same philosophy. Equations (2.51) and 
(2.52) remain the same. The only modification is the augmented 
space unit. For example/ in the 3CP.A the augmented space unit 
is the 24-gon (Figure 2,2c) (Mookerjeo/ 1973) and 
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^1 


L 2 


^2 



This has four independent components/ but note that both 
(2.51) and (2.52) for this case also yield four equations: 
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(2. S3) 
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Lbc^BC) = K 


® i 

BG ® 


^ - BM 


km 


^k 


km —mC 


As the cluster size increases/ the augmented space 
graph itself becomes large and the augmented space unit which 



rciioriual i.'.-'.oo the cluster also becomes large. For clusters of 
size 1 - 1 , thio unit is a (N x 2 ) - gon. However, the number of 
independent members in ^ still remains relatively small and the 
G can be calculated by the recursion method with great facility 
oven on quite large augmented space units. This fortunate 
marriag>.i of the above method with the recursion technique reduces 
the calculation of quite large CCPAs to tractable problem. 

The self -energy £ , hence the CCPA, is cluster diagonal, 
a property essential to preserve Herglotzicity (Butler, 1973). 

The effective Hamiltcaiian is invariant with respect to cluster 
translation. However the exact effective Hamiltonian is transla- 
tionally invariant. In contrast, in the CCPA's the self -energy 
at the cluster-centre and cluster-periphery are indeed different. 
Unless it is possible to construct clusters in which all the 
sites are equivalent (as in the 2CPA) this distinction is unavo- 
idable and inherent in any cluster approximation (Tsxikada, 1972), 
'The equations (2,51) and (2,52) form the basis of the 
self -cons latent calculation of the self -energies and 

Thu jthilosuphy behind the procedure is the standard one employed 
in all self -consistent coherent iiKSdia calG\ilations . The aim is 
to find an ”effGctive" Hamiltonian 

4*ff ^ ^ ^ (2.54) 

such that [w I - 'hhe configuration averaged resol- 

vent. Ite ^££ is found by determining the way in which the 
disorder 'renormaliaes ' diagonal term to ^(w) and off-diagonal 
term to ^ , (w ) , by immersing the cluster conf igioration in a 



modJ-Uiri and ensuring that there is no extra scattering . The 

procedure / involving a self-consistent iterative method/ is only 

useful if the iterations converge. We may steirt with the MCPA 

/ 

results as the zorootii approximation/ i.e, S°(w) = 2 „^(w) I 

^ ~ ^ these in (2.51) "'to generate and 

These in turn may bu now used in (2.52) to generate £.q(w) and. 

1 

(w ) , The procedure is then iterated until there is conver- 
gence i.e. self-consistency is achievad. 

Now we will apply this scheme on one -dimensional (Section 
2,12) and on three-dimensional systems (Chapter III), Before 
thiS/ we will discuss spectral functions and density of states 
evalxjiated by our scheme. 


2,11 Spectral Functions and Density of States 

The spectral function may be immediately deduced from 
(2.54). In K-space 


-.e£f 


N 


Z 

ij 


i.(1t.r/-k* .r / ) 


I + I, (W , ^ (k) = I X £ (w, X) 

■**o *^X V 


Wc' “O a. * 

The effective Hamiltonian is diagonal in the K-indices 


and therefor® the spectral function is 

(kO'" 


A(k/ w) 


1 r 2 - 
Im L w I 


io “ 4' 


(2.55) 


The inverse being taken over the 3 x 3 space of the longitxadinal 
and transverse branches. The spectral functions are peaked at 
the natural frequencies of the alloy and the width at half 
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maximum (FWHM) is proportional co the inverse life time of the 
st'.*to with that .jnorgy- For the CCPA/ the widths are branch 
dependv^nt .and also dependent on the direction of k. This is in 
direct contrast to the MCPA, where is diagonal and = 

0_(k) is real or tne imaginary part of h'“ (k) is branch and k- 

inde pendent. Thus the spectral functions allow us to construct 
the dispersion curves for the alloy and the accompanying widths. 
Maxima in widths indicate tsxistence of resonance y while double 
peaked spectral functions indicate strong hybridization* We 
shall extract those infoarmation for the alloys under considera- 
tion to compare v;ith neutron scattering experiments (Chapter III). 

The configuration averaged density of states is given by 
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F , (JiE) 
t ^1 


^ “* B' 

^2 = i^ol 


r V B‘M‘ M' M‘B 
+ I I K . G ^ 
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^ jk ^1 
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(2.57) 


On th»_ r .-iioruuilizud ic'-ttico/ now wc havo two states i^£-' 
p^r situ. I,--’t us indicate thorn by £ and 
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:£f 


I ' Lo Pi + d Pi ’ + ? 'df;j£ 


■P df,, Jf. "■ -2 "■ ’ 


(2.58) 


< if . j if • > rnay now be calculated from (2.58) by recursion 

starting from 1 o} = 1 o > and performing at each fraquenoy. 

o 

2*12 Results on Linea r Chains 

Here we present the results of the linear chain for a 

■ x_i. TArnr-'k Wo chooso tho mass and 

preliminary comparison with ecrlier 

Off-diagonal pur-.motors in such a way that we may compare our 

^sults with those of Mylos ^d ^ow (1979a,b) and Diehl et al. 

. • 1 K v-*-* +-h' h it is not nocBssaxy to 

(1979)# It must bo oniphasised hurc. 

i. T 1 ■ snecial values/ but any arbitrary 
limit our ciiOJ-c.- to tn^^e specia 

choice is possible- 

Figure 2.3 shows the results for vary small mass and 

Off-diagonal disordors. The results are very similar to all 

. Tn Piaure 2.4 we ahow results for mass disorder 

earlier works . In Figure 

raiir- T-Psults wlth thosc of Mylcs and Dow 
only/ and we may compare our result 

g ^ notices the following points: most 

(Figure 7)* Ohe immediat- y 

Of the fine psah structure is alxmxst identical for clusters of 
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the same size. Significant differences only occur near the band 
edges. The band width in the self-consistent field is greater 
than that given by the ICPA environment. Moreover, the sharp 
cut-off introduced by the artificial iCPAenvironment is absent, 
as well as the sharp unphysical psak there introduced by the 
squeezing in of the states into the MCPA band. It is obvious 
from the trends that larger clusters will reproduce the detailed 
peak structure more accurately, and the position and structure 
near the band edges also improve with the size of the cluster 
taken. 

Figures 2.5 and 2.6 are to be compared with Diehl et al. 
(1979, Figures 3c and b) respectively. The underlying histogram 
is that calculated "exactly" by Dean (1961) and Payton and 
Visscher (1968) for chains up to 100,000 atoms. Again, the 
method presented here is a significant improvement on the trave- 
lling pseudofermion generalization. The band edges are much 
better reproduced, as well as the structure near the band edges, 
especially in those cases where there is an impurity band with 
consider iblc structure . Agreement with the "exact histogram 
results is not bad, and it is obvious that the two will converge i 

as the size of the cluster, is increased. * 

We have put forward a formalism for calculation of ; 

averaged density of states with effect of random clustering. The 
formalism retains the analytic properties of the resolvent, whxch ; 
is crucial to a satisfactory density of states calculation and 
does so for all ranges of disorder and ani^rgy. Moreover, it : 

provides a tractable method of obtaining the self-consistent ' 
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medium/ which ia a natural -and unambiguous extension of the single 
site CF ideas. The formalism is <Mjuiklly applicable without much 
calculationnl complexity to three-dimensional systems. We are 
now in .v position to apply the formalism to actual physical 
situations: for example/ to FCC transition metal alloys in which 
much interest has been shown by experimentalists (Tsmoda et aj./1979; 
Bosi et -il/ lyfiO) duo to largo mass ana force constant disorders 
(Chapiter ill). 



Chapter III 


Application of CCPA to FCC Structured 
Transition Metal Alloys 


3.1 Introduction 

The preliminary results of calculation for linear chains 
discussed in the previous chapter were encouraging and provided 
us with sufficient theoretical basis to apply the method of 
self-consistent cluster CPA to more realistic situation of 
transition metal alloys e.g. Ni-Pd^ Ni-Cr and Ni-Pt. The choice 
is not arbitrary. In fact, these alloys form a continuous 
series of FCC solid solutions. The experimental results of 
vibrational properties of these alloys are available to compare 
with theoretical results. The mass change is dominating in 
Ni-Pd alloys? force constant change is large in Ni-Cr alloys 
due to large varicition in effective charges, while, Ni-Pt alloys 
have largt' variations both in mass and force constants (Table 1). 
This sequence makes o\jr choice for a systematic study of alloys 
with all types of disorders. The general augmented space formu- 
lation discussed in Chapter II, is developed, in particular for 

for nearest neighbour interactions only. Thus the FCC struc- 
tured alloys have another advantage of being the applicants, as 
we can neglect otlier than 1 nearest neighbour force constants 
in them without losing much structural features of physical 
properties (Chapter I } . We note here that the inclusion of 
second and third near -neighbour force constants is, in practice, 
not difficult, as we are xxsing the recvirsion technique to 
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Table 1 


fjeru,fral Properties of fcc Structured Metals (Ni^ Pt/ 
Pd and Cr) 



1 

; Ni 

-J. 

“T 

Pd 

1 

i Pt 

; Cr 

1 . Atomic 1 J umjeer 

..28 

46 

78 

24 

2. Atoi.dc mass (aii;u) 

58.71 

106.4 

195.09 

51.996 

3, Free atuin c^',s. 

3d*^4s^ 

4d^^ 

5d 6s 

0^5, 1 

3d -y.s 

4, Lattice const .ant (fcc) 
at room temp, (A° ) 

3.524 

3.8904 

3.924 

3.68 (fcc) 
(2.89) (bcc) 

5. Elastic Cs^nstcints at 

296®K (1012 dyn/cm) 

^11 

2.46 

2.270 

3.467 

3.5 

*^'12 

1.50 

1.759 

2.507 

0.678 

^44 

1.22 

0.717 

0.765 

1.01 

6. n-n force constants 
used (dyn/cm) 


18474 

22469 

27987 

37483 


19887 

22185 

31139 

17453 


-1396 

-3726 

-6663 

-13229 


TablB 2 s 

Recursion- 

-Coefficients 

for FCC 

(Ni) Lattice 


N 

■ '^IM 

r t 

1 ^2 ; 

; ‘’n i 

1 1 - ... - 1 

N 

1 1 

1 1 

i "n ■ 

L- 


1 

3.998 

4,011 

7 

4.031 

4.257 

2 

4,453 

3,548 

8 

4.176 

4.145 

3 

3.663 

3.816 

9 

3.998 

4,029 

4 

4.113 

4.579 

10 

4.085 

4.283 

5 

4.192 

3.993 

11 

4.081 

4.028 

6 

4.011 

4.092 

12 

4,050 

4.176 
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C''ilcuiatc th*.- varlouij Green's functions. The developiient of 
tho u:«j of ra.-il cpcico symmetries (Kelly ^ 1980) makes such inclu- 
sion comj(.>utcitionally possible. 

Neutron inelastic scattering experiments have been 
oxt^jnsively carried out for the vibrational properties of metals 
and alloys. The phonon line shapes/ spectral functions/ scatt- 
ering cross-sections / width functions are among the experimental 
quantities which rirC/ in general, studied. These physical quan- 
tities depends either upon G(w) or £(k' 'w)* our aim will be 
to evaluate configurationally averaged G(w) and hence G(k/ w) in 
MCPA (Eq. 2.16) and 2CPA (Eq, 2,47). Later we compare theore- 
tical dispersion relations/ width functions/ spectral functions 
and phonon density of states with experimental ones. 

Tbs theoretical calculations for G(w) in MCPA require 
the G^_(w) of pure metal, the mass ratio parameter and the 
impurity concentration (see Eq. 2.22). In 2CPA calculations we 
require aptirt from the quantities mentioned above/ off-diagonal 
elements of G^j of pure metal and force constant matrices for 

tho host and impurity. Since our present 2CPA formalism is 
applicable In the first nearest-neighbour approximation; so we 
do all calculations of pure metals and MCPA in the same appro- 
ximation in order to be consistent. Before considering any 
particular alloy we shall first discuss the force constant matrices 
for the host and impurities with the G (w) of piare metal. 

3.2.1 Force Constant Matrices 


The experimental results for the dispersion relations 
of 'gang of four' metals - Ni/ Pd/ Pt and Cr - are available 



72 


(Biry^iK.'aucjt al, 1964; Dutton et al, 1972; Muller and Brockhouse, 
1971; Peldmnrin, 1970). Those results are used/ in general, to 
fit for 3—3 U'-jiglibour force constant models to reproduce other 
physical gui.iatitiv,s e.g. density of states, and Debye-Waller 
factor t.'tc. Wo have se'„-n that for these metals except Cr which 
is hCC in metallic form at normal temperature, the nearest 
neighbour force constant is very large in Born Von Karman models. So 
we choose nearest neighbour force constants model. In this model 
theoretical 'JlSi^orsion curves are given by equation (1.45). Using 
experimental points in Eq. (1.45) we find the force constants as 
a best fit to them- The values of best fit force constant 
matrices for these metals arc given in Table 1 while the theore- 
tical dispersion curves with experimental points are given in 
Figure 1-6 .and Figures 3. 1-3. 3 for Ni, Pt, Pd and Cr metals 
respectively. It must be remembered that actual Cr is bcC. 

However Cr in dilution with Ni forms a FCC solid solution. 

Those results uiro of a hypothetical Cr FCC solid. So the compa- 
rison of the theoretical dispersion curves with the experimental 
points is not justified. 

3*2.2 Lattice Green ' s Functions 

The ]<now ledge of G^(w) of a pure metal is best obtainable 
from the Eq. (1.44) which require k-space calculations. Very 
accurate Ic— space methods exist for the calculation of tho Green s 
functions, for example, the ray integration technxque (Chen, 

1977) which was further extended by Choudhry (1981) to deal with 
off-diagonal Green's functions. Green's functions with vacancies 






Longitudinal Component 



C , Reduced Wave Vector 

FIG.3-3 Dispersion Curves for FCC Structured Cromium 
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(as ruc|,iu rcid in oiir iiiathodology ) may also be derived from these 
by isoiatud iiupurity techniques* But this requires high preci- 
sion coniputatiuriri}. t^ccurdcy* Moreover the frequent use of 

G(w)'s in s..li consistent equations (2.5l) and (2.52) makes it 
urifeasible in practice . 

Instead there is a very different approach for the 
calculation of Green’s functions — ’Recursion' technique in 
real space (ilaydock et al/ 1972; 1975) (discussed in Appendix A). 
For comparison we draw the phonon density of states for Ni 
metal with varying number of continued fraction coefficients. 
Figure 3.4 (coefficients are given in Table II). The 

k-space calculated density of states are given in Figure 3.5. 

The only difference between recursion density of states and 
k-space density of states is that near Von-Hove singularities 
recursion density of states gives oscillatory oehaviour- This 
nature can lx. smoothoned out if the asymptotic part of the 
continixd frsiction coefficients is extrapolated in an efficient 
manner. RuCuntiy thw asymptotic summation of continued fraction 
has disousse'd by JKoily (private oomniunication and SSP, 

Vol. 3S). To gut the continued fraction coefficients much 
faster onu should exploit the real space symmetries (Kelly, 1980) 
but we so far have not been able to achieve this in particular 
for phonon case (though for electronic case these are being 
exploited in our group by Joshi (private communication)). Our 
interest/ at present, is to include off-diagonal disorder effects 
.in abovo alloys. We can leave these problems as futirre exten- 
sion of this work. 



(arbitrary scale) 



1 


1 


I 



79 


H.iving known theao input parameters we are, now in 
position to study thvjse alloys in MCP and 2CP approximations. 
I’irst wo will discuss in detail Ni-Pd alloys, having large mass 
disorder . In tliu next section wo discuss Ni— Cr alloys where 
force constant disorder is aominant. Finally we wild discuss 
Ni~Pt alJ.oys whore both disorders are important. In the last 
section wo will conclude the results, 

3-3 Mi~Pd Alloy 

We choose the Ni^^Pd^^ alloy for the application of 
2CPA de»volop-*d in Chapter II because this alloy has already been 
studied theoretically (Kamitakahara and Taylor, 1974) in MCPA 
and experimentally (Kamitakahara and Brockhouse, 1974) by 
inelastic noutron scattering. The propearties of Ni and Pd have 
been jumnuirized in Table 1. Ni ana Pd form solid solutions of 
all concentrations and there are no indications of long-range 
order .Jown to 0°C. The phase diagram of that 

the liquddus anu solidus touch at 45% Pd (Hansen, 1958). The 
dispersion curves of Ni and Pd are dominated by large, almost 
central nearest neighbour forces (Birgenea.u et al, ,1964; Muller 
and Brockhouse, 1968; 1971). Force constants in Pd are about 
15% hJLgtor than Ni. 

In Figure 3.6 we display the spectral functions obtained 
from 2CPA formulation in the highest symmetry directions ( £00 
[CCOJ, teesj )- The first thing to note is that the spectral 
functions are (in contrast to Lorentzian shapes) often 
asymmetric near the , resonances . The^ asymmetries 
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occur near the reL,onance frequencies, m other words, the 
shape or a inode with a frequency slightly lower or higher than 
that ui a resonance tends to have a second peak or wide tail 
ovt.r the. resonance legion. Experimentally, for some neutron 
groups corresponding to transverse phonons with frequencies just 
belcx-/ tho lowar rosonunce , definite asyirimetries to the right 
wo*ro observed. .Such asymmetries are evident for fcczl T and 
(iJCO ! branches (Kamitakahara and Brockiiouse, 1974). In most 
casos the energy resolution is typically 0.5 THz and coianting 
statistics . ire probably inadequate to detect other types of 
.u.yiiaa' t.rios. Secondly, the spectral functions have a pronounced 
branch dependence in contrast with the MCPA predictions. 

In Figure 3.7 we display the nattiral disorder-induced 
widths (a consoquonco of spectral function of Figure 3,6) as a 
function of frequencies with the experimentally observed points 
from incoherent nautron scattering data by Kamitakahara and 
HrocKiiouse, Expurlmo'ntally, there is a evidence of resonance, 
which is strongest in tnc T branen and at 3 THz, strong 

in tlie C^CO, branch also at 3 THz. It is weaker but still 
evident in the T 2 and Coo?] T branclnes at slightly greater 

than 4 THz. In the I’eeol T 2 branch there is evidence of a 
Second resonance at around 7 Triz. There is some indication of 
two rv~ sonnnci_-h in the [c?cJ L branch but within experimental 
error tho width could increase monotonically from origin to the 
zone boundary. In the [OO?:] L branch, there is no indxcation 
at all of resonances and the amplitude of the width is smaller 

than for the other branches. 

Theoretically MCPA calculations predicts a branch 
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independent width shape because of diagonality and K-indep- 
endence of the MCPA self-energy The 2CPA calculations 
introduce both non -diagonality and k-dependence of ^ (2.55) 
and the predicted v/idth shapes are strikingly similar to experi- 
mental results (Figaro 3.7) except in T and [pOS^ L 

directions. In our calculations we do not got large width in 

T direction as experimentally observed. Though for ?< 0.25 
there is very go<.:)d matching- In [do^J L direction we still get 
two resonance widths while experimentally no resonance exists. 

But Bosi et al. (1978) got a high resonance width in L 

direction for the Pd-Pt alloy which has almost same mass ratio 
(1.836) as Ni-Pd alloy and hence, to some extent, supports onr 
prediction. 

The deviation from MCPA results is partic\ilarly notice- 
able in the concentrated alloy (45% of Pd). In these concentra- 
ted ^regimes, there is a high probability of clusters of Pd 
randomly occurring in the* alloy, these clusters have an effect of 
introducing k-dependenco and off-diagonal J^'s (in fact the off- 
diagonal has tha some matrix structior* ris the force constant, 
dictated by the symmetries of the lattice (Chapter I ) ) . For 
this alloy, such an effect also arises from disorder in the 
force constants which cannot bo taken into account in MCPA. The 
2CPA is the lowest order approximation in which it is possible 
to consider them consistently. It is difficult to guess which 
of these two effc^cts dominate in the width function although the 
analysis of Ni-Cr alloy with dominant off-diagonal disorder 
indicates a predominant effect of the latter. 
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In Figure 3.8 we plot the MCPA and 2CPA dispersion 
curves. The 2CPA further strengthens the evidence of a second 
resonance in l^5 03 T 2 branch at higher rrequencies. The kin^ in 
the transverse branches are also pronounced. Next v^e plot (in 
Figure 3.9) the u-u correlation function p(w) given by 

p'(w) = - ^ iTig « Tr dm G (1, 1; w)) » (3,1) 

in MCPA and 2CFA. In a pure crystal this is the true 

density of states (see Chapter I). But it is the natural f\anc~ 
tion to extract from incoherent neutron scattering cross-section 
results. Specific features appear in 2CPA, The enhancement of 
the peak aroxind 7,8 THz is the main feature. We may mderstand 
this as dxoe to impurity - cluster bands pulled out of the mean 
crystal density. Secondary structures appear around 3.5 THz 
where the resonance occurs. The effect of force constant dis- 
order is not transparently apparent in Figure 3.9 except perhaps 
in slight lowering the band edges in reference to MCPA calcu- 
lation . 

3.4 Ni-Cr Alloy 

Theoretically it is not a good system to study because 
of BCC structure of pxira Cr metal. But when it is alloyed with 
Ni it forms a continuous series of FCC alloys upto 30% of Cr. 

In Ni-Cr system^ the components have masses which are close 
enough to eliminate any appreciable effects of mass disorder while 
at the same time the number of outer electrons (d and s electrons 
outside the A-core) changes from six to ten and hence the effect 
of force-disorder dominates. Moreover/ Nig^Cr^^ is experimentally 
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studied by Bosi et al. (1980) to see the effect of force constant 

disorder. That is why we have chosen this system for the appli- 

cation of our theory. Our calculations are on Ni-Cr model 
system in which Cr is treated as FCC metal because of the requi- 
rement of our theory. We construct a hypothetical FCC Cr in 

such a way that its BCC elastic constants are related to FCC 
force constants in nearest neignbour approximation; 


Cfi + 2 C ^2 - 

4(f^ - f^, - f^)/a 

^11 “ ^12 = 

(fj_ + 5f^, + 2f^)/a 

''44 = 

(fj_ + \, + 2f^)/a 

The values of 

'^44 taken from Leibfried and 


Breuer (1978), The calculated values are listed in Table 1 and 
the dispersion curves are given in Figure 3,3. 

Having done this, we carry out the MCPA and 2CPA calcu- 
lations and display the some spectral functions (Figure 3.10), 
dispersion curves (Figure 3.11) and the frequency distribution 
function (Figure 3.12). In this case the widths are small in 
all directions . This is expected, as clear-cut resonances are 
features of large mass disorder. There is evidence of Ic and 
branch dependence, as also asymmetry in certain directions. 

This strengthens our arguiment that force constant disorder is 
responsible for the asymmetry of line shapes and ultimately for 
their strong k- and w-dependence , The same effect can be seen 
in the dispersion curves and in the frequency distribution 
function. The dispersion curves are lowered and distorted 

2, 5 -4, 5 THz range in all directions. Ihe lowering of 


n&iArr 




bution for 




92 


dispersion curves near edges of symmetry directions is observ- 
able in the maximum limit of density of states which is reduced. 
j i M ii i i ', i .8Hiia» gJBy . Near the broadening region (5-6 THz) of spectral 
functions, frequency distribution givas a peak while at lower 
(higher) values of this region the distribution is reduced 
compared.to MCPA. 

3.5 Ni~Pt Alloys 

Being encouraged by the right trend of theoretical 
results towards the experimental results in the Ni-Pd and Ni-Cr 
alloys , where either of two disorders — - diagonal and- off-diagonal 
dominates, we now apply our foirmulation to Ni-Pt alloys where 
both disorders are predominant. The melting points of these 
alloys are nearly constant over a wide range of concentration and 
hence tliese alloys also fojrm a continuous series of FCC solid 
solutions. The extra characters of these alloys are their atomic 
ordering at 3:1 and 1:1 compositions and hence are under the 

influence of atomic short range ordering (Hansen, 1958). The 
nearest neighbour force constant of Pt is roughly 50% larger than 

that of Ni, Tsunoda et al, (1979) have studied this system 
thoroughly covering the full range of concentration by inelastic 
neutron scattering. We, therefore, discuss this in detail for 
different concentrations (c = 0,05, 0.30, and 0.50) successively. 

3,5,1 NigsPtQs 

In Figxxce 3.13 we plot the dispersion cuirves in the 
[poO b aud [poCl directions for MCPA and 2CPA. The insets 
are enlarged for both components which explain the effect of 




Wav« Vector 

ersion curve* tor Ni^s PIqs 











0-0 0-2 04 0^ 0-8 0*0 0-25 

< , Reduced Wave Vector 
F1G.3“15 Dispersion Curves for NiggPtoB 


06 


f 

t 


t 

j 



environmental (disorder. The portion A of MCPA is stretched to hig- 
her valiaes of C in 2CPA, while portion B of MCPA is pushed up to 
higher frequencies in 2CPA. This pushi-ig and pulling is the 
consequence of hybridization of ideal ^phonon lines with the 
resonance mode which is possible only if S. has strong k- and w- 
dependence (Dederichs and Seller, 1380). The ultimate effect is 
the increment in the energy gap with the resonance frequency. 

The anamolies corresponding to the resonance modes of Pt occurs 
at nearly the same frequency of 3,05 THz in MCPA for both compo- 
nents, while it is near 3.32 THz in 2 CPA calculations. The 
experimentally observed frequency mode is 3.3 THz almost equal ; 

to our prediction. At this frequency both components have 
double peak structxire. Some of the spectral functions (line 
shapes) for the (]00?3 L branch are given in Figure 3.14 after the 
calculation by MCPA and 2CPA. For other two directions L^^ol | 

and dispersion curves are given in Figure 3.15. In MCPA 1 

i 

the resonance frequency appears to be insensitive to the phonon 
propagation direction. In 2CPA we also obtain sarre nature except | 
in h direction (where resonance mou.e occurs at about 3-27 j 

THz). Since at this low concentration one cannot, from experi- j 

j 

f 

mental results, claim the dependence or independent of k- vector j 

! 

in the resonance mods, so we shall later use our formulation for j 
higher concentrations to see the directional dependence effects, j 

I 

in Figure 3.16 we display the MCPA and 2CPA frequency distribu- j 

i 

tion fxanctions along with the partial density of states xn the | 

inset calculated by Tsunoda et al. by inelastic neutron scattering 

! 

f 

experiment • The experimental partial density of states are . 
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increased in frequency range 4-6 THs. Similar effect can be 
seen in frequency distribution function in the same region for 
2CPA. Moreover, second peak is also increased and near the 
resonance frequency there is a kink in 2CPA. As far as the 
maximum frequency is concerned tnere is not much change, in 2CPA 
calculations except a slight lov^ering of the band width. 

3.5.2 Ni7QPt30 

The dispersion curves along the [00 sj/ and 

directions are displayed in Fig\ire 3.17 for MCPA and 2CPA. The 
resonant behaviour observed for c = 0.05 is strongly enhanced in 
this system. For the [OO T branch the phonon spectral f xanct- 
ions (line shapes) becones well separated double peaks in the 
region from C - 0,15 to the zone boundary in MCPA; from C= 0,30 
to the zone boiandary in 2CPA; while from C = 0.40 to the zone 
boundary experimentally. For C>’0,5, both the MCPA and 2CPA 
indicate that the 00 C T branch increases with ? to the zone 
boundary, wher-as the experimental branch seems to decrease in 
this region. This leads to a difference in interpreting the 
position of the resonance mode. Tsunoda et al. interpret this 
as the mean between the maximum energy at C 0.5 on either 
side of which [OO Q T branch decreases, and the lowest energy 
of the upper toOCj T branch. This gives a value of 3,5 THz. 

But in the MCPA and 2CFA we interpret this as the itean between 
the maximxjm of the lower T branch at the zone boundary and the 
minimum of the upper T branch. It is 3.75 THz in MCPA while 
4.1 IHz in 2CPA. However, if we take the maximxim of the lower T 
branch at around S - 0.45 beyond which the peak value of spectral 
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function is decreasing^ we shall get 3.6 THz a mean resonance 
value in the 2CPA in good agreement with the experimental value. 
Another feature in the 2CFi\ result is a sudden shift upward in 
the lower spectral peak aroxind ? 0.45 in the T branch and 

similar shifts in other branches ( [cco3 end L branches). The 
widths here are very large and the peak values themselves low 
so that individual interpretations become difficult to compare. 

For the Coo i' branchy MCPA calculation predicts the 
existence of a resonance mode at v v 7 THz in addition to the 
main resonance at ^ 3.75 THz but at the secondary resonance 

the calculation shows only an asymmetric line shape with a large 
width. In 2CPA calculation we obtain large hybridization at 
4.1 THz ( ? = 1.5-^. 5 IHz) and smaller secondary splitting 
near Z 0»65 and v 6,8 THz while observed splits are near 
C 0.35, V 4.0 THz and ^ 0.6, v 6.0 -mz. Near secondary 

resonance 2CPA shows an asymmetric spectral functions of large 
width but with a dip resulting in the splitting. In fact, as 
discussed earlier in Ni-Pd case, such a secondary resonance is 
due to more scattering induced by impurity at higher frequencies 
causing resonance modes to decay into the host bands with the 
stxrength being proportional to the host phonon density of states 
and hence increasing the line widths. 

The dispersion curve of C^^o3 T^ has the almost same 
nature as [oo T except that resonance is at about 3.65 THz in 
2CPA. The [^?o3 L component has one resonance mode only at 
^ = (1.0 - 2.5) mm and ^ 13,9 IHz in 2CPA. Ihe absence of 
second resonances in this branch and ^ xinderstandable 

if we consider the maximum limit of impurity density of states 





which is approximately 6 THz. The introduction of large concen- 
tration of impurities will affect the Ni-Ni force constant much and 

♦ 

hence the high frequency range of pure density of states (6-8 
Ttiz) with the result than a second resonance appears around 6 THz, 
But at higher > 0.4} in these branches the higher frequencies of 
host atom dominate not being affected very much by impurity 
atoms. We do get second resonance in T 2 branch which is 

not observed in experiments. In T branch MCPA and 

2CPA Show double peak structure in almost entire region one 
below '^ 3.5 THz and one above it. The resonance frequency is 
'^3.4 THz in MCPA/ '^3.5 THz in 2CPA while 3.5 THz is experiment- 
ally observed. In ^ branch 2CPA resonance width decreases 

( '^'0.4 THz) while in MCPA width is 0.5 THz). 

In Figure 3,18 we display the frequency distribution 
function for MCPA and 2CPA with the experimental partial density 
of states in the inset. The experimental observation shows the 

t 

sudden increment in density of states at about 4.0 THz (a 
primary resonance mode region) and pushing in thtj band limit at , 
higher frequencies. The same narvire we get by 2 CP^ calculation. | 
There is sharp enhancement of frequency distribution at about I 

4.0 THz and at about 7.6 THz and band is also reduced. In a | 

partial dGnsxty o£ states for the host atom calculation/ the j 

[ 

impurity peaks near the band edges will he removed. The remainder j 

1 

of the frequency distribution function in the 2CPA is in csxcellent ; 
agreement with the experimental curve shown in the inset. 

The 2CP/» calculations show the branch dependence of 
resonance mode in better way than by MCPA. The trend being m ! 
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the right direction is encouraging. Next we study the 50-50 
• * 

alloy , 


3*5.3 NisQPtsQ 

We display the dispersion curves along the symmetry 
directions in Pigixre 3.19^ after MCFA and 2CPA calculations. 

MCPA results do not show direction oi' branch dependence. In 
[ oof} L, t??0l [OOSiT and branches resonance occurs 

at about 4.15 tHz in MCPA calculations. It is about 4.4 THz in 
other directions e.g, L, t???! T and [c? ol L. In ^ 

branch the width is largest compared to other branches. It is 
similar to the nature observed in Wi^^Pd^^. In 2CPA calculations 
there are few extra feafures: the secondary resonance is enhanced 
in [00?^L and (5^0j T 2 directions, the L is more distor- 
ted* there is prominent kink in TjL T direction 

the branch and direction dependence is easily visible , in 
■[00 5 J T direction resonance occiors at about 3.95 THz (if we 
follow the procedure mentioned in the previous section ) , at ? = | 

0,6, The experimental value is 3,5 THz. The primary resonance ; 

! 

in [pO Cj b direction occurs at about 4,9 THz cind the secondary ! 

one is at about 6.8 THz. Experiment shows only one resonance of | 
large width at about 5,2 THz. Although the experirresntal inteS""* i 
pretation is difficult in this region, because of large width ; 

(Tsunoaa et al, 1979). ; 

Now we compare the freguency distribution function, £qk 1 

i 

MCPA and 2CPA. The nature found in 2CPA is same as in the previoxis 
alloy supporting our views. Experimental partial density of I 

states are shown inside the inset. 2CPA results follow the trend . 
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of experimental obseirvatlorB near the resonance mode and at the 

band edge. In 2CPA we get one extra small delta function like 

peak near the band edge. The origin of this peak is not very 

that 

clear to us. It seems to bo^he some frequency of host atom is not 

affected by the presence of impurity (an off-diagonal disorder 

effect)/ or may be due to cluster effects. Again partial density 

of states will not show these cluster structure of impurities 

in 

and the rest of the frequency distribution is^excellent agreement 
with the exporiment_(inset)/far superior to the MCPA calculation. 

3.6 Conclusion 

In conclusion it must be said that despite the 2CPA 
being first approximation with off-diagonal disorder, its appli- 
cation to FCC structured transition metal alloys is found to be 
qviite satisfactory qualitatively and quantitatively as well. In 
all the three systems considered above, most of the observable 
quantities are better explainable by 2CPA as compared to the 
existing MCPA predictions , The k- and w-dependence of spectral 
functions and width functions is an achievement of the 2CPA. 

Since nearest neighbours are large in number in FCC lattice and | 

I 

environmental disorder is properly accounted for in 2CPA, the j 

effect of larger clustering is not expected to be important in I 
these systems, I 


I 



Chapter IV 


Coiicludlnq ' Remarks 

A seif -consistent cluster coherent potential approxi- 
mation (CCPA) is developed for the study of disordered alloys of 
materials very different in masses and force constants. This is 
a significant development over existing MCPA (ICPA). The formu- 
lation, being general, is applicable to any kind of elementary 
excitations. Out of CCPA's, 2CPA is fully developed for the 
study of vibrational properties of substitutionally disordered 
binary alloys. 

To sum up, the following are the achievements of this 

work: 

(1) A self-consistent CCPA based on augrrtejnted space formalism 
of Mookerjee, has been developed to take into accoxint the 
effect of large clusters, 

(2) The translational invariance condition is properly incor- 
porated which seems to be difficult to include in conven- 
tional CPA theories (Mostoller and Kaplan, 1979; Grunewald, 
1976). 

(3) The analytic properties (Herglozicity ) of the Green's 
functions are maintained in the CCPA theory at all energies, 
concentrations and degrees of disorder. 

(4) Out of CCPA's, 2CPA theory is applied to PCC structured 
transition metal alloys in the systematic way - with 
diagonal disorder only; off-diagonal disorder only and 
with both disorders being predominant. 



However, there are a few refinements still essential 
over the calculations presented here . They are s 

(1) The real-space methods for the calculation of different 

2 

Green's functions (G. .(w )) need further refinements. 

Since we go down up to 12 steps in recursion method, there 
are sorre spurious oscillations near the von-Hove singu- 
larities. Those oscillations can be avoided by going down 
more steps in the recursion. This can be facilitated by 
the use of real space symmetries (Chapter I) which will 
reduce the working time enormously. Haydock (Solid State 
Physics, Vol, 35) has shown that the entire basis 1 u^^ > 
of tigljt -binding functions is unnecessary, for the recur- 
sion method, only those cembinations which have the 
symmetry of the lattice and the starting state lo) 
contribute . This drastically reduces the "effective" basis 
and thus CPU time on the computer. 

Another approach is to take a finite cluster and use 
the 'super smooth' programme of the Cambridge group to 
smoothen out the resulting histograms. 

(2) We have considered nearest -neighbour force constant inter- 
action obtained from Born-von Karman models. Instead we 
can use other near neighbours also for the exact input 
information. This can also be achieved by reciarsion and 
some modifications in our formalism. Such alteration will 
provide a general formulation for the applications to 
other lattices e.g, SC, BCC and diamond. 
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(3) We can solve the configurationally averaged Green's 
function (2.46) alternatively. Instead of following the 
above self -consistc-nt scheme ^ we can exploit graphical 
method with recursion directly. On these lines some 
initial work has been carried out by Hersovici and Fibich 
( 1980a, b), in which they have considered only diagonal 
disorder. With all disorders we are trying to exploit 
-this scheme, 

(4 ) Though we have concentrated mainly on the determination of 
spectral functions, width functions, dispersion curves and 
frequency distributions, it is also possible to extend 
this formalism for the calculation of vibrational density 
of states, cluster density of states, etc. used in the 
calculation of scattering cross-sections, elastic constants 
and thermal conductivity of alloys with little modifications 

(5) In this theory, we coiild include Markovian type short-range 
order (Srivastava, 1982). 

(6) The effect of bigger clusters can be included by performing 
a recursion calculation in the full augmented space. 

So far, we have confined oiorsolves only to crystalline 
systems with compositional disorder.. Amorphous alloys, where 
there is a departure from crystalinity, but there is an under- 
lying distorted network on which the atesn sit (Meek, 1976; 

Yamamoto et al, 1980) may also be considered by the use of 
'recursion* with the tractable size of a cluster. These re»>„ 
finements we leave as possible future extensions of this work. 
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In addition, as an extra piece of investigation, we 
have studied the screening density around a static charge in an 
electron liquid using a dielectric function proposed by Tripathi 
and Mandal (1C77) (See Appendix B). The non-linear behaviour of 
this function is due to the approximate averaging scheme of 
double -time retarded Green’s function in their Eq. (2.14). This 
calculation is an important check of tne accuracy of their 
approximation in evaluating the dielectric function. 
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Appendix A 


Recursion Method 


The 'recursion* method of Haydock et al. (1972; 1975) 
is developed for finding density of states as well as other 
matrix elements of the resolvents and wavefunctions. Its moti~ 
vation was that since many observable properties of amorphous 
and disordered solids depend mostly on local structure^ it ought 
to be possible to formulate a theory in which local effects are 
explicitly separated from bulk averages. It is closely related 
to Lanczos method (Wilkinson ^ 1965). The philosophy of tiiis 
technique is as follows. For a given Hamiltonian matrix in 
a suitable Hilbert space the recixrsion method generates a 
unitary transformation such that 

U H (A.l) 


where H^^^ is tridiagonal i.e. 


H * 

— TD-imn 


^m+l 

Hr- 


n = m 
n = m -f 1 
n = m - 1 
otherwise 


(A. 2) 


For this^ new basis functions are constructed iteratively so 
as to guarantee that each new member interact with its preceding 
and following members. The first fxonction in this new basis is 
chosen to be the one whose local density is sought. Here/ we 
ncM f discuss its practical use. 



The representation of the Hamiltonian H for the vibra- 
tional system given by oq. (1.1) is suitable for the use of 
recursion method. Our interest is in finding the one-phonon 
spectral density given by (1.46). Rewriting it > 

p(w) = - ^ I Im (w^ + is ); ( 6 = O"^) 

TT a 


where 

oo ,2, , , 2 _ 1, 

^aa / oj(w j[“2^ |a/0> 


(A. 3) 


with 


D = M 


2 m-1/2 


M and ^ are the mass-tensor and force constant matrix, respec- 
tively. la , o > is the component of the vibrational state of 
the atom o. p (w) is normalised (in modes per atom) as 


w 

.max 

i P (w ) dw = 3 


where w^. is the maximxam frequency. 

rncix 

. oo , 2 V . 

For constructxng G (w ) as a continued fraction we 
start with .i. given clement > = 1 a, o > and an arbitrary sequ- 
ence of numbers define a sequence of elenents 

{| u^> } by 




D i u 


h-l 


a - 

n-l 


tu 

' n-l 


- b' 


n-l 


l^n-2" 


= 0 


(A. 4 ) 


The sequence termixiates at n if null vector (b^ = 0) 

othexTwise we generate an infinite sequence of elements. The 
numbers real if H(D) is Hermitian, If the sequence 

of elements { | } is such that <u^ I ^ n ^ m then 

equation (A. 4) immediately implies 



< u 


n 


D 


n 


<u 


n 


<u, 


*n 


n--l 


D 


u > 
n 


<u 


n 


u > 
n 


n 


n 


<u. 


‘n-1 


1 > 
n-“l 


<u 


n-l 


u - > 
' n-1 


(A*5) 


If we take {ju^>} as our basis of representation the Hamiltonian 
D becomes tridigonal with the down the diagonal and 

down the off-diagonal positions • 

If D is a non-Hermitian matrix (as is often the case 
with effective Hamiltonian arising out of any coherent potential 
approximation) then, alternatively we choose an arbitrary set 
of complex numbers a^/ b^, as follows 

bo |Ui> = (D-aoHUo> - =0"“ 

for n > 0 


n 


u > 
n-1 


(2 - a^) 


u > 
n 


c 


n 


u 

n- 


1 


> 


Using Equation (A*6), it is easy to see that 


n 


‘^u D u > , 
n ‘ n 


n 


< u ^ D u > 
n-l ^ ‘ n 


(A. 7) 


“b “ ''"n I ‘2 - (D - a„) |u_^> - |c„| ; 


if we assuxne <u 1 u > = . 

n ' m nin 

The off-diagonal elements of the Green's operator 

11 ' 2 

, (w ) may also bo calculated in the same way. For example, 
if we wish to calculate 

G°g (w^) = < 1 (w^ I - D)"^ 

then we must start with 


Uo> 


1 

2 


[a0> + 1&«1> 

' «* ' 3 ? 



which givec 


0 

1 (W^ 

1 " D)”'^ 

1 

' 0 

But 

00 

^aa 

00 

= 0^3 and 

pOl 


for configurational 
homogeneous in space. 


1 

2 


[G 


oo 

aa 


+ G 


11 

$3 


+ G 


ol 

a3 


+ 




~ general, for pure lattices and 

Green's functions in case the disorder is 
So 


Gap ) 


= <u„ ! (w I 


2)" 




- G 


OO 
a a 


(w^) 


(A. 8) 


Hera the first term of the right hand side of above 
equation can be calculated by continued fraction expansion with 

+ 1^ .,1> and Gqq is the term which has already 

been obtained by the same procedure with [n^> = |a^0>. Hence 

ol 2 11 ’ 2 

we can calculate G^g (w ) . Similarly other (w ) are also 

calculated. 

The reason why this method is called 'recursion ‘ msthod 
now becomes clear. Starting from a Hamiltonian and an element 
ju^> with = 0 we may generate with the help of equation (A*4), 
a^ / b^ and thus i > . Now we may determine in turn a^ , b^ and 
thus j U 2 > and so on recursively. As a calculation technique 
this method is very economical. At each stage we need to 
store and |u^ and the coefficients calculated upto that 

stage. All extraneous information is left out. ; 

2 ^ ' 

Now, let M be the determinant of the matrix (w I - : 

and (n = 0, 1, 3, ...) be the determinant of the matrices ! 

derived from (w^ I - D^) after rows and columns from 0 to n : 

have been eliminated. Then if we wish to obtain (w ) = 

i 

we may use the Cauchy expansion for the determinant and write | 
down a recursion relation for as | 



(A. 9) 


= (w ~ \+2 

oo 2 

•^aa ' becomes a continueci fraction expansion from the repeated 
application of equation (A, 9) 


(w^ - a^) » b?M. 

O O 11 


(w^ “ ■" 


ai -b^ 


"2 “ ^3 


(A. 10) 


, r 2 -,2 , 2 . \ ; 

« - ^'n-l - ‘'n I 

I 

f 

By the recursion method we can calculate iteratively [ 

i 

(a / b ) from n to some N. The effect of all further iteratronsi 
' n n^ o i 

can be absorbed in if Ca„y bj) converge fast (which in I 

r L , ^ n n 

^ ”^2 2 2 ^ 

general/ happens). Thc;n ^ ^ ^ will give j 


CO C3D 


an analytical expression for the remainder resolvent 


= G (w ) 

^N+1 


2 >.2 

W d., . ^ ““ -D*, 


■^+1 N+2 

^ ~ ^+2 ” ^N +3 


(A,ll) 


w^ - a - (w^) 

CO CD CD 


Which yields 



A- 6 


G (w ) 
00 ‘ 


1 2 


(w^ - a ) - - a - 4b^ ] 


(A. 12) 


We reject tho second solution of the quadratic as the 
) must V‘ .nish as w" -»■ oo in order that the resolvent pres- 
erve its analytic character (of the form there. It is 

clear that ) is real for real outside the range 

“ 2 ] ) and + 2 ! )/ s^ that the continued 

fraction will be real. Inside the same interval G (w^) has a 

ou 

branch cut and we obtain a band of allowed energies. Since 
cannot be negative (a consequence sum rule), its lower limit 


is 0, hence a 


2b ^ holds in the phonon density of states 


calculation by recursion. If wo taiow w homehow (here we used 

ni6oc 

K-space maximum frequency limit in Born -von Karman model). 

We can fix exact value of a , h and hence G (w"^ ) . 

oo ciD ou 

Now we can obtain g'^'^(w^) from (A. 10) and hence phonon 
density of states 


(w) = 


^o "^1 


"1 ~ ^2 


, (A. 13) 


The tridiagonalisation of D thus reduces the original 
system to a fictitious linear chain of atoms with local vibra- 
tional displacements with diagonal energies a^ and n-n inter- 
actions b^ (Solid State Physics, Vol. 35). 
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Department of Physics, Indian Institute of Technology, Kanpur^ ^ 
Dielectric Function, Static Screening, and Non-Linear Response 
By 

R. P. SINGH and M. YUSSOUFF^^ 


Tnpathy and Mandal /l/ proposed a new dielectric function e (k,C£> ) solving 
the equations of motion of the double-time retarded Green^ s functions by de- 
coupling the higher order Green’ s functions G(k,a) ) with the conservation of 
frequency moments up to infinite order. The resulting G(k, a? ) satisfies the 
compressibility sum rule and yields the asymptotic value of G(k) equal to 1/3. 
Unlike previous theories /2, 3/ here G(k) has a sharp peak around k = 2k^. 

The new features of e (k, co ) given by Tnpathy and Mandal arise from 
the approximations made in evaluating it. These authors proposed a decoupling 
scheme which amounts to expressing the higher order Green’ s function in 
terms of the lower order Green’ s function along with a frequency and wave- 

number dependent coefficient . (k,Ci:) ) in their equation (2.14). But their 

^ 1 ’ 1 

evaluation of this coefficient is rather arbitrary because they fall back on the 

same equation involving the higher order Green’ s function to evaluate A-j ^ (k,a3). 

1 

Then they evaluate the former approximately by trying to retain terms in 

A-*- . (k, oo ) that are linear in v(k). The effect of such an unclear approxi- 

^1 ’ ^1 

maiion may be expected to show up in the static limit which reflects itself in 
the screening density around an impurity of unit charge. The following calcu- 
lation shows that this is the case. 

In the linear response theory, the screening density at a distance r from 
the static impurity of unit charge situated at the origin is given by / 4/ 


&Q (r) 


2x 


oc 

/ 


dk 


ksin(k r) 


e (k, 0) 


). 


(1) 


where r is in units of k'^^ and k in units of k^. The static dielectric function 
according to Tripathy and Mandal is given by 


1 - 


Qo(k) 


ilkTO) =1 +Q (k)[l -G(k)T 


( 2 ) 
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where 


2 ar 


Q (k) = ^Fik) 


(3a) 

with 

F(k) =1 + i(l - 4 k )ln| k. 2 1 ' 


(3b) 

substituting ( 2 ) in ( 1 ) one gets 



6 D(r) 1 ksin(kr) 

.3 ~ 0 ^^ X ^ ^ ' 

k » 2 sc 0 

Q(j(k) 

(4) 

Q^(k)Ci - G(k)T 


The numerical values of G(k) I 



0 Oi 08 U Iff Z.0 LB 32 ?i 



given by Tripathy and Mandal are 
used to evaluate the integral ( 4 ) using^ 
Simpson’s rule. Very good conver- ^ 
gence was obtained with the interval ^ 
0. 01 to the upper limit lOO. The re- j 
suits are compared with those of 
Pathak /5/, Singwi and Tosi / 6 /, anc^/ 
Toigo and Woodruff / 7/ for Tg = 3 / 

in Fig. 1. / 

It is clear that the present resultj 
are very different from the earlier 
ones. The amplitudes of oscillations 
are larger by almost a factor of two ! 
and the detuiy is slow. But the initial ; 
value of 6 g (r) is close to that ob- ■ 
tained in earlier work. All these can I 
be explained in terms of the behaviou^ 
of G(k) as k varies. Since the asymp- 
totic limit of G(k) is comparable to 
the Hartree-Fock value the behaviour of 

3 

Fig, 1 . Screening density 6 g(r)Af 
versus rk, for r = 3; a) for small 
values of ^ rk- and b) for large values 
of rk-. (1) Singwi and Tosi / 6 /, 

(2) Pathak /5/, (3) Toigo and 
Woodruff /7/, (4) present results 
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5§(r) for small r is close to earlier formulations. Also G(k) of Tripathy and 
Mandal is similar to that of others in the range of O^k^l which governs the 
jecay at large r values. The difference in the range IS k^ 2 with a notable 
eature that a peak exists around k = 2k^ is responsible for the large amplitude 
iscillations and decay at moderate values of r. Thus the value of G(k) around 
= 2 k^ is debatable and the above calculations might serve as an experimental 
heck of the accuracy of Tripathy and MandaP s approximation . 

One of us (R. P.S. ) is indebted to CSIR, New Delhi, for a Junior Research 
'ellowship while we were working on this problem. 
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